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I.  INTRODUCTION  AND  PURPOSE 


The  goal  of  the  research  effort  reported  here  has  been  to  derive 
a  class  of  stochastic  attrition  models  from  probabilistic  assumptions 
on  the  behavior  of  individual  combatants  and  on  the  interactions  among 
them.  Our  interest  is  in  stochastic  analogs  of  a  family  of  determin¬ 
istic  attrition  models  commonly  called  ’’Lanchester  attrition  models”, 
after  their  originator  F.  W.  Lanchester.  These  deterministic  models 
are  discussed  in  Section  II. 

There  are  several  reasons  for  undertaking  such  an  effort. 

Previous  research  on  stochastic  Lanchester  models,  as  surveyed  in 
Section  III,  has  been  concerned  mostly  with  certain  computational 
problems  and  in  general  imposes  by  hypothesis  the  form  of  the  attri¬ 
tion  process,  rather  than  deriving  that  form  from  more  elementary 
assumptions.  In  some  cases,  therefore,  our  derivations  lead  to  known 
and  studied  processes;  the  point  is  that  instead  of  arbitrarily  imposed 
processes  we  deal  with  the  consequences  of  elementary  and  physically 
meaningful  hypotheses.  In  the  terminology  of  economics  we  employ  a 
’’micro”  rather  than  a  "macro”  approach,  stating  assumptions  about 
individual  combatants  rather  than  about  the  overall  form  of  the 
attrition  process. 

There  should  be  a  general  preference  for  stochastic  rather  than 
deterministic  attrition  models.  A  stochastic  model  is  more  general, 
more  flexible,  more  realistic,  and  better  founded,  and  always  provides, 
through  expectations  of  its  outputs,  scalar  characterizations  of  the 
system  being  modeled.  Deterministic  attrition  models  of  Lanchester 
type  are  represented  by  differential  equations  and  require  allowing 
noninteger  numbers  of  combatants,  while  the  stochastic  models  we 
present  here  have  only  integer  (but  vector-valued)  states. 
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The  reasons  for  wishing  to  derive  models  from  elementary  assumptions 
are  also  several.  First,  one  wants  to  know  if  there  exists  a  set  of 
assumptions  from  which  a  known  process  can  be  derived  and,  if  so,  what 
physical  situations  are  consistent  with  the  assumptions.  Understanding 
of  the  model  and  its  possible  applicability  are  enhanced  when  assump¬ 
tions  are  explicitly  stated.  Different  models  of  combat  can  be  com¬ 
pared  in  a  reasonable  manner  on  the  basis  of  underlying  assumptions 
(as  well  as  by  their  relationship  to  historical  data)  rather  than  on 
the  untenable  basis  of  outcomes,  implications,  and  heuristic  judgments. 

Once  sets  of  assumptions  are  in  hand,  new  models  may  be  created 
by  generalizing  or  weakening  certain  assumptions.  For  example,  we 
have  found  that  a  certain  stochastic  Lanchester  model  frequently  used 
for  describing  combat  between  heterogeneous  forces  is  not,  based  on 
underlying  assumptions,  the  appropriate  generalization  of  the  corres¬ 
ponding  homogeneous  model.  The  effect  of  weakening  unrealistic  and 
untenable  assumptions  can  be  explored  in  a  sensible  way  only  if  it 
is  realized  what  those  assumptions  are. 

Finally,  once  underlying  assumptions  are  found  for  a  family  of 
related  processes,  one  may  seek  general  structural  characteristics, 
unifying  taxonomies,  and  general  computational  approaches  such  as 
we  present  in  Section  IV. 

The  final  section  of  this  paper  is  a  compendium  of  the  processes 
derived  so  far,  giving  for  each  process  one  family  of  assumptions 
from  which  it  can  be  derived  and  a  probabilistic  characterization  of 
the  process.  It  may  be  that  certain  of  these  processes  can  be 
derived  from  alternative  sets  of  assumptions,  but  we  believe  that 
our  families  of  assumptions  are  essentially  unique. 

Probabilistic  technicalities  and  proofs  of  our  results  appear  in 
the  Appendix. 

The  stochastic  attrition  processes  discussed  here  are  all  time- 
dependent  dynamic  models  with  a  continuous  time  parameter.  In  Karr 
(1972a,  1972b,  1973,  1974)  similar  derivations  are  given  for  a 

class  of  static  and  discrete  time  attrition  models. 
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II.  DETERMINISTIC  LANCHESTER  ATTRITION  THEORY 

Consider  a  combat  between  two  homogeneous  forces,  Blue  and  Red, 
and  denote  by  b(t)  and  r(t)  the  numbers  of  Blue  and  Red  survivors  at 
time  t  after  the  combat  is  initiated.  The  British  engineer  F.  W. 
Lanchester  (1916)  suggested  that  it  is  the  nature  of  modern  warfare 
that  the  instantaneous  casualty  rate  on  each  side  be  proportional  to 
the  current  strength  of  the  opposing  side.  Lanchester  thus  proposed 
the  now  famous  model 


(1) 


b'(t)  =  -  c^t) 
r'(t)  =  -  c2b(t) 


where  c^,  c ^  are  positive  and  not  necessarily  equal.  In  order  that 

(1)  make  sense,  the  functions  b  and  r  must  be  allowed  to  assume 
arbitrary  nonnegative  values. 

The  solution  of  (1)  subject  to  the  initial  conditions 

b(0)  =  bQ 
r( 0 )  =  rQ 

is  given  by 

b(t)  =  b^  cosh  Xt  -  sinh  Xt 

(2) 

r(t)  =  r^  cosh  Xt  -  a  b^  sinh  Xt 

where 

\  =  (cic2^ 
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and 


a  =  (Cl/c2^ 


X  is  a  measure  of  the  intensity  of  the  engagement  and  a  of  the 
killing  effectiveness  (per  unit  time)  of  one  Red  combatant  relative 
to  that  of  one  Blue. 

The  functions  b  and  r  defined  by  (2)  are  of  interest  only 
until  the  time  t  =  infjt:  b(t)  =  0  or  r(t)  =  0} at  which  one  side  or 
the  other  is  annihilated,  t  is  infinite  if  and  only  if 


in  which  case 


r(t)  =  lim  b(t)  =  0 


lim 

t-400 


t— k® 


otherwise  one  side  is  annihilated  at  a  finite  time  and  the  other  has 
a  positive  surviving  strength. 

It  follows  from  (2)  that 


2r  2  ...2-,  ,2  , /t,2 

a  [rQ  -  r(t)  ]  =  bQ  -  b(t) 


(3) 


for  all  t  .  In  view  of  (3)  the  system  (1)  of  differential  equations 
is  called  Lanchester’s  "square  law"  of  attrition. 

Lanchester  also  proposed  the  so-called  "linear  law"  in  which  each 
side’s  casualty  rate  at  any  time  is  proportional  to  the  product  of  its 
strength  and  the  strength  of  the  opposition.  In  differential  form 
this  model  is  given  by 


b'(t)  =  -  k1b(t)r(t) 


(4) 


r'(t)  =  -  k2r(t)b(t) 
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where  k-^,  k^  are  positive  constants.  These  are  not  the  constants  of 
(1)  and,  indeed,  must  have  different  units.  The  exact  solution  to 
(4)  does  not  interest  us  (but  note  that  it  is  nonnegative,  removing 
one  objection  to  (2)).  The  condition  analogous  to  (3)  is 

Tq  -  r(t)  =  constant  x  (b^  -  b(t)) 

from  which  the  term  "linear  law”  originates. 

Of  considerable  interest  is  the  distinction  between  ’’square 
law”  and  ’’linear  law”  combat.  Historically  [see,  for  example, 

Bonder  (1970,  p.  160),  Deitchman  (1962,  p.  818),  Dolansky  (1964, 
p.  345),  and  Hall  (1971,  p.  8)]  the  belief  has  been  that  the  square 
law  describes  combat  situations  in  which  individual  opponents  are 
identified  and  engaged  one-by-one,  a  situation  commonly  referred 
to  as  ’’point  fire”  combat.  On  the  other  hand,  the  linear  law  has 
been  thought  to  describe  combat  processes  in  which  weapons  (such  as 
artillery)  fire  only  at  an  area  in  which  opponents  are  located 
(’’area  fire”).  The  research  presented  here  leads  to  a  number  of 
relevant  conclusions,  whose  overall  effect  in  the  preceding 
distinction  is  somewhat  ambiguous.  Indeed,  perhaps  the  best 
distinction  is  the  following.  ’’Square  law"  and  "linear  law”  are 
terms  which  distinguish  two  types  of  engagement  initiation.  In 
the  former  one  side  initiates  engagements  (i.e.,  fires  shots)  at 
a  mean  rate  which  is  proportional  to  its  own  numerical  strength 
and  independent  of  the  numerical  strength  of  the  opposition,  while 
in  the  latter  each  side  initiates  engagements  at  a  mean  rate  pro¬ 
portional  to  the  product  of  its  numerical  strength  and  that  of  the 
opposition.  On  the  other  hand  the  "area  fire”  -  "point  fire" 
distinction  seems  to  involve  mostly  what  can  happen  once  a  shot 
is  fired.  Area  fire  processes  allow  multiple  kills  with  one  shot 
while  point  fire  processes  allow  at  most  one  kill  by  a  single  shot. 
This  dichotomy  is  not  unambiguous;  for  example  certain  processes 
which  seem  physically  to  represent  the  "firing  at  an  occupied  area" 
aspect  usually  associated  with  area  fire  may  in  fact  have  at  most 
one  kill  per  shot,  due  to  dispersion  of  one  side  or  the  other. 
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Seen  in  this  way  the  square  law- linear  law  and  point  fire-area 
fire  are  not  competing  distinctions  but,  rather,  are  complementary, 
so  that  in  fact  four  classes  of  processes  are  defined  by  this  classi¬ 
fication  scheme. 


Let  us  now  consider  the  case  of  heterogeneous  forces,  where 
each  side  consists  of  two  or  more  distinct  types  of  combatants. 
Suppose  there  are  M  types  of  Blue  weapons  (combatants)  and  N 
types  of  Red  weapons  and  denote  by  b^(t)  and  r^(t)  the  number  of 
Blue  type  i  and  Red  type  j  weapons,  respectively  surviving  at 
time  t  .  Let  b(t)  =  (b^t),  ...,  bM(t))  and  r(t)  =  (r^t),  ..., 
r^(t))  be  the  vectors  of  Blue  and  Red  surviving  forces  at  time  t, 
respectively. 


Since  there  is  no  known  set  of  precise  and  unambigous  assumptions 
leading  by  means  of  a  rigorous  mathematical  derivation  to  (1)  (Weiss 
(19S7)  is  typical  of  the  ambiguity  and  vagueness  of  previous  sets 
of  assumptions  and  the  lack  of  rigor  of  previous  derivations)  it  is 
not  clear  what  should  be  the  appropriate  generalization  to  the 
heterogeneous  case.  The  usual  heterogeneous  version  of  (1)  is 
obtained  purely  formally,  by  replacing  b,  r,  c^  and  c 2  there  by  the 
vector  b,  the  vector  r  and  matrices  c^  and  c^  to  obtain  the  model 


(5) 


N 

b((t)  =  -  E  c  (i, 
j=l 


j )rj ( t) 


M 

r'(t)  =  -  E  c2(j,  i)b.(t) 
J  i=l 


i  =  1,  .  .  . ,  M 

D  d  5  •  •  •  J 


which  is  often  abbreviated  in  matrix  form  as 
b'  =  -  c  r 
r'  =  -  c2b  , 

where  c^  =  [c^(i,  j )]  and  c2  =  [c2(j,  i)]  are  M  x  N  and  N  x  M  non¬ 
negative  matrices,  respectively. 
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Evidently  there  exist  neither  a  closed-form  solution  to  (5)  nor 
any  means  of  handling  negative  components.  The  termination  rule  so 
easily  described  for  the  homogeneous  square  law  model  (1)  does  not 
carry  over.  Another  of  our  discoveries  is  that  the  stochastic  attri¬ 
tion  process  whose  analog  is  (5)  is  not,  based  on  sets  of  assumptions, 
the  appropriate  generalization  to  the  heterogeneous  case  of  the  sto¬ 
chastic  process  analogous  to  (1).  Hence  the  common  practice  of 
calling  (5)  ’’heterogeneous  Lanchester-square  combat”  is  unjustified. 

As  part  of  our  research  we  have  derived  the  proper  stochastic  and 
deterministic  heterogeneous  square  law  models,  which  are  presented 
in  Section  V. 

Relatively  little  has  been  said  about  heterogeneous  versions  of  the 
linear  law  model  (4)  although  it  has  been  suggested  by  D.  Howes  (per¬ 
sonal  communication  to  L.  B.  Anderson)  that  the  appropriate  form  is 

N 

b'(t)  =~bi(t)  E  k1(i,  j)r.(t) 
j=l  J 


(6) 

M 

r'(t)  =t  .(t)  E  k2(j,  i)bi(t) 
J  3  i=l 


where  k  ,  k^  are  nonnegative  M  x  N  and  N  x  M  matrices,  respectively. 
What  the  motivation  for  (6)  is,  other  than  (4)  and  intuition,  is 
not  clear,  but  it  does  turn  out  that  the  stochastic  version  of  (  6) 
is  the  proper  generalization,  based  on  assumptions,  of  the  stochastic 
version  of  (4). 

The  symmetry  of  (1)  and  (4)  is  not  necessary.  One  might  assume 
instead  that 


(7) 


b'(t)  =  -  cr(t) 
r'( t)  =  -  kb( t)r( t)  , 
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for  positive  constants  k  and  c  ,  which  is  the  so-called  "mixed- 
law”  of  Lanchester  combat.  The  mixed  law  was  proposed  by  Deitchman 
(1962)  as  a  model  of  guerrilla  ambushes,  with  Blue  representing  the 
ambushers  and  Red  the  ambushed. 
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III.  SOME  ASPECTS  OF  PREVIOUS  WORK  ON 
STOCHASTIC  LANCHESTER  PROCESSES 

We  do  not  intend  here  to  survey  in  detail  the  large  collection  of 
previous  studies  of  stochastic  models  purported  to  be  the  proper  analogs 
of  the  deterministic  systems  (1),  (4),  (5),  (6)  and  (7),  except  by 
indicating  the  philosophical  and  mathematical  differences  between  those 
approaches  and  our  approach  here.  The  two  most  extensive  published 
surveys  are  Dolansk#  (1964)  and  Hall  (1971);  Springhall  (1968)  also 
contains  a  survey  and  a  rather  extensive  bibliography.  Included  in 
the  list  of  references  are  a  number  of  papers  of  interest  in  the 
historical  development  of  Lanchester  theories  of  combat;  in  particu¬ 
lar,  the  papers  of  Brown  (1955),  Isbell  and  Marlow  (1956)  and  Weiss 
(1957)  contain  the  germs  of  the  theory.  Sets  of  assumptions  in  these 
early  works  are  imprecise  and  ambiguous  and  derivations  are,  for  the 
most  part,  incomplete  or  nonexistent.  One  purpose  of  this  research 
is  to  provide  consistent  and  unambiguous  sets  of  assumptions  from 
which  certain  stochastic  attrition  models  can  be  rigorously  derived. 

In  doing  so  we  have  unified  and  extended  the  theory  of  Lanchester 
combat  models . 

In  order  to  understand  the  attrition  processes  presented  in 
Section  V,  some  background  concerning  a  certain  class  of  Markov 
processes,  namely  regular  step  processes,  is  required;  we  shall  also 
make  use  of  the  following  discussion  in  the  remainder  of  this  section. 
Blumenthal  and  Getoor  (1968),  Freedman  (1971),  and  Karlin  (1968)  are 
principal  references. 

Let  E  be  a  countable  set.  A  Markov  matrix  on  E  is  a  mapping 
P  :  E  x  E  -*  [0,  1]  with  the  property  that 
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E  P(i,  j)  =  1 
jeE 


for  all  i  e  E. 

Given  a  Markov  matrix  P  on  E  such  that  P(i,  i)  =  0  for  all  i 
and  a  function  X  :  E  -  [0,  «)  there  exists  a  Markov  process  (Xt)t  q 
with  state  space  E  satisfying  the  following  intuitive  description. 
If  \(i)  =  0,  then  once  the  process  enters  state  i  it  remains  there 
forever  after,  while  if  X(i)  >  0  the  process,  upon  entering  state  i  , 
remains  there  an  exponentially  distributed  time  with  mean  1/X(i) 
independent  of  the  past  history  of  the  process,  whereupon  it  jumps 
to  another  state  in  E  according  to  the  probability  distribution 
P(i,  *)j  independent  of  the  length  of  its  sojourn  in  state  i.  (X^_ ) 
is  called  the  regular  step  process  with  jump  function  \  and 
transition  kernel  P  . 

The  family  (P  of  Markov  matrices  on  E  defined  by 

PtU,  :>  =  P{xt  =  j|x0  =  i| 

is  called  the  transition  function  of  (X^).  The  matrix- valued  mapping 
t  —  p  can  be  shown  to  be  differentiable  and  the  matrix 


has  the  property  that  P^  =  QP^  for  all  t  .  Q,  is  called  the 
infinitesimal  generator  of  the  regular  step  process  (X^)  and  is 
given  by 


(8) 


Q(i,  j)  = 


-  Hi) 

X(i)P(i,  j) 


if  j  =  i 
if  j  1  i  . 


Q  has  the  interpretation  of  specifying  the  "inf initesimalff  or 
"differential"  behavior  of  the  process  (X^)  because  for  j  /  i, 
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Q(i,  j  )  is  the  "infinitesimal  rate”  at  which  the  process  (Xt) 
moves  from  state  i  to  state  j  in  the  sense  that 

PtXt+h  =  j'Xt  =  il  =  ^i>  3)  h  +  °(h) 

as  h  -*  0.  Here  lim  o(h)/h  =  0.  The  resemblance  to  a  system  of 
hiO 

differential  equations  is  evident,  so  a  given  stochastic  attrition 
process  is  called  an  analog  or  version  of  one  of  the  deterministic 
models  (1),  (4),  (5),  (6),  or  (7)  provided  its  infinitesimal 
generator  sufficiently  resembles  the  appropriate  system  of  differ¬ 
ential  equations.  Further  details  are  in  the  Appendix. 

For  example,  the  second  stochastic  attrition  process  presented 
in  Section  IV  has  infinitesimal  generator  Q  given  by 

Q((i,  j);  (i,  j  -  D)  =  icl 

(9)  Q((i,  j);  (i,  j))  =  -  (icx  +  jc2) 

Q((i,  j);  (i  -  1,  3))  =  jc2 

where  c^,  c^  are  positive  constants  derived  from  quantities  given  in 

the  appropriate  family  of  assumptions.  The  first  equation  in  (  9  ) 
says,  in  the  differential  interpretation  of  Q ,  that  when  Blue  and 
Red  strengths  are  i  and  j  ,  respectively,  Red  casualties  (that 
is  the  transition  from  j  Red  survivors  to  j  -  1  Red  survivors)  are 
occurring  at  infinitesimal  rate  ic^  Hence  there  is  justification 
for  calling  the  process  whose  infinitesimal  generator  is  the  Q  of 
(9),  a  stochastic  homogeneous  square  law  attrition  process,  because 
of  the  clear  resemblance  between  (1)  and  (9).  This  justification, 
incidentally,  is  far  from  new,  dating  at  least  to  Snow  (1948). 

All  previous  work  on  stochastic  Lanchester-type  processes  begins 
essentially  at  (9)  by  imposing  as  a  hypothesis  the  form  of  the 
generator  of  the  process  to  be  studied.  Such  studies  have  generally 
been  concerned  with  computing  quantities  of  interest  such  as 
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(1)  expected  numbers  of  survivors  at  each  fixed  time; 

(2)  distribution  and  expectation  of  the  time  required 
to  reach  certain  subsets  of  the  state  space  (such 

as  the  set{(i,  j):  i  =  0  or  j  =  Of  of  absorbing  states 
which  is  entered  when  one  side  or  the  other  is  anni¬ 
hilated); 

(3)  the  probability  that  Blue  wins  the  engagement  by 
exterminating  Red. 

Our  concern  has  been  directed  at  a  more  fundamental  problem.  It 
is  elementary  to  show,  for  example,  that  there  exists  a  regular  step 
process  whose  infinitesimal  generator  is  the  Q  given  in  (9),  but 
no  one  has  previously  presented  a  complete  and  unambiguous  set  of 
assumptions  on  the  firing  behavior  and  interaction  (or  lack  thereof) 
among  combatants  which  entail  an  attrition  process  whose  infinitesimal 
generator  is  the  Q  of  (9).  It  is  this  lack  of  basic  and  physically 
meaningful  sets  of  underlying  assumptions  that  this  research  attempts 
to  alleviate. 

Hence  in  some  cases  the  processes  we  derive  are  known,  but  not 
always.  In  all  cases,  it  is  the  derivation  and  the  more  primitive 
level  of  the  underlying  assumptions  which  are  new.  Sometimes  one 
set  of  assumptions  leads  the  way  to  a  new  process  (in  particular 
this  is  the  case  for  the  stochastic  versions  of  the  heterogeneous 
square  law  (5))  and  these  new  processes  are  also  discussed  in 
Section  V. 
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IV.  STRUCTURE  OF  THE  FAMILY  OF  PROCESSES  DERIVED 


We  discuss  in  this  section  several  unifying  structures  and  taxon¬ 
omies  which  can  be  applied  to  the  family  of  stochastic  attrition 
processes  presented  in  the  next  section,  from  which  the  reader  hope¬ 
fully  can  obtain  both  overview  and  insight. 

The  first  taxonomy  classifies  the  stochastic  attrition  processes 
presented  here  in  terms  of  three  criteria: 

(1)  Multiple  kill  (area  fire)  or  single  kill  (mainly  point  fire) 

(2)  Square  law  engagement  initiation  or  linear  law  engagement 
initiation 

(3)  Homogeneous  or  heterogeneous  force  compositions 

This  appears  in  Table  1  below.  The  classification  as  to  type  of  en¬ 
gagement  initiation  is  based  on  analogy  between  infinitesimal  genera¬ 
tors  of  the  process  and  the  various  systems  of  differential  equations 
appearing  in  Section  I. 

In  each  instance  the  homogeneous  model  is  a  special  case  of  the 
heterogeneous  model  and  each  heterogeneous  process  correctly  reduces 
to  the  corresponding  homogeneous  process  when  each  side  consists  of 
only  one  type  of  combatant. 

The  second  taxonomy  is  based  on  the  families  of  assumptions  under¬ 
lying  the  processes,  which  are  presented  in  Section  V.  The  taxonomy 
appears  in  Table  2,  which  gives  for  each  of  the  processes  the  form  of 
the  basic  assumptions  from  which  it  is  derived  (except  assumption  (3) 
which  holds  for  all  of  the  processes).  A  typical  family  consists  of 

(1)  An  assumption  that  either 

(a)  Times  between  shots  fired  by  a  surviving  weapon  are 
independent  and  identically  exponentially  distributed 
with  some  mean,  that  when  a  shot  is  to  occur  exactly 
one  opponent  is  detected,  attacked,  killed  or  not,  and 
lost  from  contact,  all  instantaneously;  or 
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Table  1.  TAXONOMY  OF  THE  PROCESSES 


Force  Composition 


Engagement 

Initiation 

Homogeneous 

Heterogeneous 

I.  Multiple  kill  (area  fire) 

Square  law 

1)  Process  Al,  a  new  process 
with  multiple  kills  and 
square  law  engagement 
initiation 

II.  Single  kill  (m 

ainly  point  fire) 

Square  law 


1)  Process  SI,  whose  genera¬ 
tor  is  analogous  to  (1) 


1)  Process  S2,  whose  genera¬ 
tor  is  analogous  to  (5) 


2)  Process  S3a,  a  new  process 
obtained  by  simple  exten¬ 
sion  of  the  assumptions  of 
Process  SI 


3)  Process  S3,  a  new  family 
of  processes  incorpora¬ 
ting  fire  allocation,  of 
which  S3a  is  a  special 
case 


Linear  law 


1)  Process  LI,  whose  genera¬ 
tor  is  analogous  to  (4) 

2)  Process  L2,  a  new  process 
generalizing  Process  LI  by 
the  inclusion  of  engage¬ 
ments  of  positive  duration 


1)  Process  L3,  a  new  process 
obtained  by  extension  of 
the  assumptions  of  Process 
LI  to  the  heterogeneous 
case,  with  generator 
analogous  to  (6) 


Mixed  law 


1)  Process  Mia,  a  special 

case  of  Ml,  with  generator 
analogous  to  (7) 


1)  Process  Ml,  a  new  process 
of  which  a  process  with 
generator  analogous  to 
(7)  is  a  special  case. 
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Table  2.  ASSUMPTIONS  OF  THE  PROCESSES 


Process 

Assumptions 

A1  Homogeneous 
Square 

la)  Mean  firing  rate  is  specified. 

2)  Shot  kills  binomially  distributed  number  of 
opponents. 

SI  Homogeneous 
Square 

la)  Mean  firing  rate  is  specified. 

2)  Shot  kills  exactly  one  opponent  with  probability 
p,  none  with  probability  1  -  p. 

S2  Heterogeneous 
Square 

la)  One  firing  process  for  each  opposing  weapon  type 
with  rate  dependent  on  target  and  shooter; 
all  such  processes  occur  simultaneously  and 
independently . 

2)  Kill  probabilities  depend  on  target  and  attacker. 

S3  Heterogeneous 
Square 

la)  Each  shooting  weapon  fires  shots  at  mean  rate 
dependent  only  on  its  type.  Allocation  of  fire 
over  opposing  weapon  types  is  by  prescribed 
probability  distributions. 

2)  Kill  probabilities  depend  on  target  and  attacker. 

S3a  Heterogeneous 
Square 

la)  Same  as  S3  except  that  fire  allocation  is  by 
uniform  distributions  (special  case). 

LI  Homogeneous 
Linear 

lb)  Mean  time  to  make  one-on-one  detection  specified. 
2)  Kill  probability  specified. 

L2  Homogeneous 
Linear 

lb)  As  in  LI. 

2)  Engagement  is  one-on-one,  lasts  for  exponential¬ 
ly  distributed  duration,  ends  with  death  of  one, 
the  other,  or  neither  combatant  with  specified 
probabilities . 

L3  Heterogeneous 
Linear 

lb)  Mean  detection  time  depends  on  target  and 
searcher. 

2)  Kill  probability  depends  on  attacker  and  target. 

Ml  Heterogeneous 
Mixed 

1)  Each  side  possesses  two  weapon  types,  one  of  which 
behaves  according  to  la)  and  the  other  according 
to  lb ) . 

2)  Kill  probability  depends  on  target  and  attacker. 

Mia  Homogeneous 
Mixed 

1)  One  side  has  weapons  described  by  12);  the  single 
weapon  type  on  the  other  side  is  described  by  lb). 

2)  Kill  probability  specified. 
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(b)  The  time  required  to  detect  a  particular  opponent  is 
exponentially  distributed  with  some  mean,  different 
opponents  are  detected  independently,  and  every 
opponent  detected  is  instantaneously  attacked,  killed 
or  not,  and  lost  from  contact; 

(2)  Specification  of  necessary  conditional  probabilities  of 
kill  given  detection  and  attack; 

(3)  An  assumption  that  firing  processes  of  all  combatants  are 
mutually  independent.  Thus  each  weapon  operates  independently 
of  all  weapons  on  the  other  side  and  all  other  weapons  on  its 
own  side. 

In  heterogeneous,  mixed,  area  fire  (Al)  and  time-to-kill  (1,2) 
models,  these  assumptions  are  weakened  or  modified.  The  heterogene¬ 
ous  Process  S3  requires  additional  assumptions  concerning  allocation 
of  fire.  In  heterogeneous  processes,  mean  rates  of  fire  or  detec¬ 
tion  and  kill  probabilities  depend  in  general  on  both  target  and 
shooting  weapon  types. 

The  universal  independence  assumption  (3),  even  though  it  is 
omitted  from  Table  2,  should  not  be  overlooked.  It  states  that 
in  a  probabilistic  sense  there  is  no  interaction  among  weapons  on 
a  given  side  and  interaction  among  weapons  on  opposing  sides  only 
when  a  kill  occurs .  In  particular,  none  of  these  models  is 
capable  of  handling  synergistic  effects  sometimes  thought  to  be 
important,  except  perhaps  by  artificial  (and  possibly  unjustifiable) 
devices  such  as  modifying  the  initial  numbers  of  weapons  of  some 
types,  based  on  the  absence  or  presence  of  some  other  weapon  type 
before  applying  one  of  the  attrition  models  presented  here. 

We  have  attempted  to  keep  our  assumptions  as  free  from 
restrictive  physical  interpretation  as  possible,  in  order  to 
demonstrate  the  full  range  of  applicability  of  each  model.  For 
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example  consider  the  assumption  (1)  of  Process  SI  (see  Section  V) 
which  states  that  times  between  shots  fired  by  a  surviving  weapon 
are  independent  and  identically  exponentially  distributed.  This 
assumption  is  compatible  with  a  number  of  different  physical  reali¬ 
zations  of  combat.  One  can  envision  combatants  as  stationary  and 
firing  at  rates  dependent  only  on  their  own  nature  (this  seems  to 
be  an  ’’area  fire’’  kind  of  assumption)  or  as  pressing  forward  in 
such  a  way  as  to  maintain  a  constant  mean  rate  of  engagements  with 
the  opposition.  The  point  is  that  our  assumptions  are  not  unique 
in  terms  of  physical  situations  in  which  they  might  be  felt  to  be 
satisfied  and  we  have  endeavored  to  state  them  in  terms  which  make 
it  easy  to  verify  their  plausibility. 

The  widely  held  and  already  mentioned  belief  in  the  correspondence 
Square  Law-*-»-Point  Fire 

Linear  Law-*-*-Area  Fire 

is  misconceived.  Rather,  as  we  have  mentioned  earlier,  there  exists 
the  following  two  way  classification  of  processes 

Engagement  Initiation 
Square  Law  Linear  Law 

Multiple 

Kill  Structure  - - - - 

Single 


which  appears  to  us  to  make  good  sense. 


We  may  look  also  at  the  common  properties  of  the  family  of 
stochastic  attrition  processes  we  have  derived.  Among  these 
properties  are  the  following: 

1)  All  are  regular  step  processes  (in  particular,  all  are 
Markovian) ; 

2)  Infinitesimal  generators  resemble  the  Lanchester  differential 
models  of  combat; 

3)  All  components  of  sample  paths  are  nonincreasing  (we  have  not 
included  provision  for  reinforcements); 

4)  All  states  are  either  transient  or  absorbing  (the  latter 
represent  extermination  of  one  side  or  the  other). 

Quantities  of  interest  one  seeks  to  compute  for  use  in 
modeling  attrition  would  include 

1)  Expected  numbers  of  survivors  (and  thus  expected  attrition) 
at  various  fixed  times  after  the  combat  begins; 

2)  The  distribution  and  expectation  of  first  entry  times  of 
various  subsets  of  the  state  space  (which  might  represent, 
for  example,  breakoff  points  or  extermination  of  one  side); 

3)  The  expected  numbers  of  survivors  at  random  times  such  as 
those  in  2)  above; 

4)  Variances  of  certain  quantities,  for  use  in  estimating  errors 
made  in  computational  implementations; 

5)  In  heterogeneous  models,  expected  attrition  caused  by  each 
type  of  opposition  weapon. 

An  advantage  of  having  a  family  of  attrition  processes  with  some 
common  characteristics  is  the  possibility  of  developing  general  com¬ 
putational  methods.  We  will  now  briefly  discuss  one  which  might  be 
used  to  compute  expected  attrition.  For  concreteness,  consider  a 
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homogeneous  process  ((B^,  R^))^q  where  B  and  R  are  the  numbers 
of  Blue  and  Red  survivors  at  time  t,  respectively,  With  initial 
conditions  of  i  Blues  and  j  Reds  the  expected  number  of  Blue 
survivors  at  time  t  is  given  by 


i  3 


(10)  E[B  |(B0,R  )  =  <i,j)]  =  Z  E  kP  ((i,j);  (k,A)) 
r  u  u  k=0  i=0  r 


where,  as  will  be  recalled, 

Pt((i,j);  (k,jQ)  =  Pj(Bt,Rt)  =  (k,£)|(B0,R0)  =  (i,  j )}  . 

But  by  known  properties  of  infinitesimal  generators  we  have 

oo  n 

PfCCi,  j);  (k,  £))  =  E  Qn((i,  j);  (k,  £)) 
r  n=0  n- 

+-Vi 

where  Q  is  the  n  power  of  the  infinitesimal  generator  matrix  Q 
of  the  process.  Therefore  (interchange  of  the  order  of  summation  is 
justified  by  boundedness  of  all  quantities)  we  can  write 

ErB.I(Bn,Rn)  =  (i,j)]  =  Z  £  E  E  k  Qn ( ( i ,  j);  (k,  £))  . 
r  U  U  n=0  n-  k=0  1=0 

OO 

In  particular,  the  infinite  sum  £  is  the  limit  as  M  -  00  of  the 

n=0 

M 

finite  sums  £ 
n=0 

Thus  for  each  M  we  have  the  approximation 


M  tn  1  0  n 

(11)  E[B  |(B0,R  )  =  (i,j)]  ~  £  £  £  k.Q((i,j);  (k ,j>)) 

u  n=0  n>  k=0  SL=0 
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When  M  =  1  this  becomes,  provided  instantaneous  multiple  kills  are 
precluded, 

(12)  E[Bt|(B0,R0)  =  (i,j)]  ~  i  -  tQ((i,  j);  (i  -  1,  j)) 

which  is  essentially  the  analogous  (deterministic)  differential 
equation.  This,  incidentally,  provides  a  justification  (albeit 
tenuous)  of  using  the  deterministic  differential  models  as  approxi¬ 
mations  to  the  true  expectations  obtained  from  stochastic  models . 

But  why  choose  M  =  1  in  (11),  once  the  whole  family  of  approxi¬ 
mations  (11)  is  available  and  justified  and  when  it  can  even  be  shown 
that  each  is  a  better  approximation  than  the  preceding  one?  The  matrix 
Q  consists  mostly  of  zeroes  Tindeed,  Q(  (i,  j);  (k,  l))  is,  in  general, 
zero  unless  (k,  a)  is  (i  -  1,  j),  (i,  j  -  1),  or  (i,  j)]  so  that  it 
would  be  feasible  to  take  M  =  3  or  4  or  5  in  (11)  to  produce  more 
accurate  approximations  to  the  true  expectations.  Moreover,  this 
scheme  of  approximation  is  valid  for  all  the  processes  we  have 
derived.  The  principal  difficulty  would  be  computation  of  powers  of 
Q  but  since  Q  consists  mostly  of  zeroes  it  might  well  be  feasible 
to  MdiagonalizeM  Q  ;  that  is,  to  find  a  representation  of  the  form 

Q.  =  A-1DA 


where  D  is  a  diagonal  matrix,  in  which  case 


.n  -l_n_. 
Q  =  A  D  A 


becomes  trivial  to  compute  and  approximations  in  ( 11)  with  M 
essentially  arbitrary  would  become  feasible.  Indeed  in  this  case, 
we  can  obtain  P  explicitly  as 


P  =  a  ^  etD  A. 


Further  research  is  desirable  on  such  matters. 
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V.  A  COMPENDIUM  OF  STOCHASTIC  LANCHESTER  ATTRITION  PROCESSES 


We  present  here  the  family  of  attrition  process  as  derived  in 
this  research  effort.  The  format  is  to  give  for  each  process  the 
set  of  assumptions  from  which  we  have  derived  it,  and  a  description 
of  its  properties  as  a  regular  step  process,  namely,  its  jump  function, 
transition  kernel,  and  infinitesimal  generator,  followed  by  inter¬ 
pretive  remarks.  We  also  give  specific  versions  of  expression  (12) 
for  each  process. 

We  first  present  some  remarks  on  the  general  structure  of  the 
families  of  assumptions  leading  to  the  processes  described  in  this 
section.  Each  family  includes  assumptions  concerning  the  following: 

1)  Structure  of  the  two  opposing  forces; 

2)  The  rate  at  which  engagements  are  initiated; 

3)  The  evolution  and  possible  outcomes  of  an  engagement  (in  all 
processes  but  L2  an  engagement  occurs  instanteously  and  can 
end  only  in  death  of  some  weapons  on  the  engaged  side  and 
loss  of  contact); 

4)  If  necessary,  allocation  of  fire; 

5)  Interaction  among  weapons  on  a  given  side; 

6)  Interaction  of  a  weapon  on  one  side  with  weapons  on  the 
opposing  side,  other  than  engagement. 

In  all  models,  5)  and  6)  are  independence  assumptions  stipulating 
no  interaction. 

As  general  notations  we  establish  B.  for  Blue  survivors  at  time  t 
and  R^  for  Red  survivors  at  time  t.  These  will  sometimes  be  vectors. 

For  each  process  we  give  the  values  of  jump  function,  transition 
kernel,  and  infinitesimal  generator  only  for  nonabsorbing  states. 


21 


The  set  of  absorbing  states  is  left  to  the  reader  to  determine  (it  is 
the  family  of  states  in  which  one  side  or  the  other  is  exterminated). 
If  x  is  absorbing  \(x)  =  0,  Q(x,  y)  =  0  for  all  y  and  P(x,  •)  is 
not  defined. 


?? 


A1  Homogeneous  Square  Law  Area  Fire  Process 


Assumptions 

1.  All  weapons  on  each  side  are  identical. 

2.  Times  between  shots  fired  by  a  surviving  Blue  weapon  are 


independent  and  identically  exponentially  distributed  with  mean  1/r^. 


3.  If  a  shot  is  fired  by  a  Blue  weapon  when  there  are  k  Reds 
surviving  the  shot  has  the  instantaneous  effect  of  killing  a  number 
of  Red  weapons  which  is  binomially  distributed  with  parameters 

(k,  Px). 

This  means  that  the  shot  has  probability  p  of  killing  each 
particular  Red  weapon  and  that  different  Red  weapons  are  killed 
independently  of  one  another. 

4.  Red  weapons  satisfy  assumptions  2  and  3  with  mean  time  between 
shots  1/r,  and  individual  target  kill  probability  p2- 

5.  Firing  processes  of  all  weapons  are  mutually  independent. 
Process  Characterization 

Under  the  preceding  set  of  assumptions  the  stochastic  survivor 
process  ((Bt,  Rt))  Q  is  a  regular  step  process  with  jump  function 
\  given  by 


\(i,  j)  =  ir1(l  -  (1  -  P-l)3)  +  Jr2(  1  -  (1  -  P2)1)  > 


transition  kernel  P  given  by 


l  =  0,  . . . ,  j  -  1 


k  =  0,  . . . ,  i  -  1 
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and  infinitesimal  generator  Q  given  by 

Q((i,  j);  (i,  jO)  =  i^f^Cl  -  px)A- '  pj_"£,  1= 0,  ...,  j-1 

QCCi,  j);  (i,  j))  =  -Cir^l  -  (1  -  PX)D)  +  jr2(l  -  (1  -  P2)1)] 

QCCi,  j);  (k,  j))  =  jr2^(l  -  p2)k  p^-k,  k=0,  i-1  • 

The  binomial  distribution  of  assumption  2  can  be  replaced  by  any 

probability  distribution  on  0,...,k  ,  so  this  process  is  in  fact  a 
whole  class  of  similar  processes.  In  particular,  the  distribution 
which  places  mass  p^  on  1  and  mass  1  -  p^  on  0,  yields  the  Process 
SI  described  below. 

This  process  is  in  many  ways  the  most  interesting  of  those  we 
have  derived,  because  of  the  number  of  intriguing  questions  it 
raises.  The  interpretation  of  its  assumptions,  especially  that 
concern  binomially  distributed  multiple  kills,  seems  to  be  unequi¬ 
vocally  that  of  area  fire.  (Indeed,  the  possibility  of  multiple 
kills  may  well  be  the  distinguishing  feature  of  area  fire  combat 
processes).  What  is  at  issue  is  whether  this  is  a  square  law 
process  or  a  linear  law  process.  Consider  its  infinitesimal 
generator  Q.  Examining  the  term  -  Q((i,j),(i,j)),  which  is  the 
rate  at  which  the  process  leaves  state  (i,j),  we  have  what  appears 
to  be  of  square  law  form  (cf.  the  corresponding  term  for  process 
SI).  But  on  the  other  hand 

Q((i,j),  (i,j  -  1))  =  irxj(l  -  P1)'i”1P1  > 

which  can  be  thought  of  as  a  representation  of  the  instantaneous 
rate  at  which  Red  combatants  are  being  destroyed  when  the  current 
state  of  the  process  is  (i,j),  can  be  argued  to  be  of  linear  law 
form  (but  not  entirely,  because  of  the  presence  of  the  factor 
(1  -  p1)3"1).  • 

One  way  to  clarify  the  situation  is  to  abandon  the  square  law- 
linear  law  dichotomy  for  attrition  entirely  and  to  simply  regard 
the  processes  we  have  derived,  in  terms  of  their  underlying 
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assumptions.  What  distinguishes  Process  A1  from  all  the  other 
processes  is  the  provision  for  multiple  kills,  so  it  might  even 
be  argued  that  this  is  the  only  true  "area  fire"  process  presented 
here.  The  whole  situation  is  complicated  by  the  fact  that  if  the 
binomial  distribution  of  assumption  (2)  is  replaced  by  the  distribu¬ 
tion  with  mass  p  on  1  and  mass  1  -  p  on  0  the  Process  SI,  which  seems 
clearly  to  be  a  square  law  process,  is  obtained. 
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SI 


Homogeneous  Square  Law  Process 


Assumptions 

1.  All  weapons  on  each  side  are  identical. 

2.  Times  between  shots  fired  by  a  surviving  Blue  weapon  are 
independent  and  identically  exponentially  distributed  with  mean  1/r  ;  . 

3.  When  a  shot  occurs  it  kills  exactly  one  Red  weapon  with 
probability  p^  and  no  Red  weapons  with  probability  1  -  p  ,  instan¬ 
taneously  and  independent  of  past  history  of  the  process. 

4.  Red  weapons  satisfy  assumptions  2  and  3  with  mean  time  between 
shots  of  l/r^  and  kill  probability  p^. 

5.  The  firing  processes  of  all  weapons  are  mutually  independent. 
Process  Characterization 

For  i  =  1,  2, let  c_^  =  r_^p_^.  Then  subject  to  Assumptions  1  to  5 
above,  ((B^,  Rt))  q  is  a  regular  step  process  with  state  space 
N  x  N,  jump  function  \  given  by 

\(i,  j)  =  ic1  +  jc2  , 


transition  kernel  P  given  by 


P((i,  j);  (i,  j  -  D) 


1C. 


1CX  +  JC 


2 


P((i,  j);  (i  -  1,  j))  =  ieiVj  c2  . 
and  infinitesimal  generator  Q.  given  by 

Q((i,  j),  (i,  j  -  1))  =  ic1 

Q((i,  j),  (i,  j))  =  -  (ic1  +  jc2) 

Q((i,  j),  (i  -  1,  j))  =  ic 2  • 
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For  this  process  the  approximation  (12)  is  given  by 
E[Bt|(B0,  Rq)  =  (i,  j)3  ~  i  -  jc2t  ; 

similarly 

E[Rtl(B0,  Rq)  =  (i,  j)3  ~  j  -  iCjt  . 

Clearly  the  larger  t  is,  the  worse  the  approximation. 

By  the  form  of  Q.  this  process  is  the  stochastic  version  of  the 
deterministic  square  law  (1).  1/c^  is  the  mean  time  between  fatal 

shots  fired  by  a  Blue  weapon,  so  c^  =  r^p^  can  be  interpreted  as  the 
mean  rate  at  which  a  Blue  weapon  kills  Red  weapons  (which  is  inde¬ 
pendent  of  the  strength  of  the  Red  force  so  long  as  the  latter  remains 
nonzero),  giving  in  physical  terms  what  the  coefficients  c^  and  in 
(1)  mean. 
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S2  Heterogeneous  Square  Law  Process 


Assumptions 

1.  There  are  M  Blue  weapon  types  and  N  Red  weapon  types. 

2.  Consider  a  Blue  type  i  weapon.  The  times  between  shots  it 
fires  at  Red  type  j  weapons  are  independent  and  identically  ex¬ 
ponentially  distributed  with  mean  1/r^j,  i). 

3.  These  N  firing  processes  occur  simultaneously  and 
independently. 

4.  A  shot  fired  by  a  Blue  type  i  weapon  at  Red  type  j  weapons  is 
directed  at  one  Red  type  j  weapon  which  is  chosen  according  to  a 
uniform  distribution  and  is  fatal  with  probability  p^Cj,  i). 

5.  Red  weapons  also  satisfy  2-4,  with  parameters  r2(i,  j), 
p^Ci,  j)  describing  the  fire  of  a  Red  type  j  weapon  against  Blue 
type  i  weapons . 

6.  Firing  processes  of  all  weapons  are  mutually  independent. 
Process  Characterization 

For  q  =  1,  2  and  appropriate  k  and  i  ,  let  Cq(k,  l)  = 
r  (k,  £)p  (k,  £).  Note  that  B  and  R  are  now  vectors  of  dimension 

H  H  L  t 

M  and  N  ,  respectively.  Then  subject  to  assumptions  1  through 
6  above,  the  vector-valued  process  ((Bt>  ?t^t>0  a  re9u-*-ar  steP 
process  with  state  space  N  x  ...  x  N  (M  +  N  times)  and,  for  states 


(x;  y)  =  (Xp 


•  •  •  9 


•  •  •  9 


y  ),  jump  function  \  given  by 


X(x;  y)  =  E  E 

i:Xj>0  j :yj>0 


(c2(i,  j)y_j  +  c1(j,  i)xi) 


transition  kernel  P  given  for  j  with  y_.  >  0  by 


P((X,  y);  (x;  y±,  . . . ,  y . -1 
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and  for  i  with  >  0  by 


R((x,y)3  ( x^ » •  • . , x^  •  •  •  > >  y ) ) 
and  infinitesimal  generator  Q  given  by 


N 

E 

j=l 


c2(i,j )y^ 
l(x,y) 


> 


M 


Q(  (x,y);  ( x  j  y1#  • .  .,yj  1 » •  •  •  >  ^  c^(j,i)x^  if  y^  >  o 


Q((x,y);(x1, 


Q((x,y);  (x,y) )=  -EE  (c2(i,j)y. 

i  :x  .>0  j  :y  .>0  -1 

i 

+  c1(j,i)xi) 


N 

x±  ~  1>**»>xm;  y))=  E  c2(i,j)y.  if  x;L  >  0. 
’  j=l 


State  (x,y)  is  absorbing  if  and  only  if  x  =  0  or  y  =  0. 

Corresponding  to  the  approximation  (12)  we  have  for  this  process 


for  i  =  1,  M  while  for  j  =  1,  . ..,  N 

M 

E[r£|(B0,  R0)  =  (x,  y)]  ~  yj  -  (  _E  CjL(j,  i)x±)t  . 

It  seems  rather  natural  to  call  this  process  "heterogeneous 
Lanchester  square"  because  of  resemblance  of  Q  to  the  system  (5) 
of  differential  equations,  and  this  process  provides  a  physical 
interpretation  of  the  matrices  in  (5).  But  the  physical  interpreta¬ 
tion  of  the  assumptions  is  quite  unappealing.  In  particular,  assump¬ 
tions  2  and  3  are  unpalatable.  They  state  that  each  weapon  carries 
out  one  firing  process  for  each  type  of  opposition  weapon  and  that 
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all  these  processes  are  evolving  concurrently  and  independently. 

On  the  contrary,  the  spirit  of  assumption  2  of  processes  A1  and 
SI  is  that  a  weapon’s  firing  behavior  is  a  function  only  of  its 
type  and  of  whether  there  are  a  positive  number  of  targets  present. 
The  next  process  we  present  is  in  that  spirit. 

It  is,  in  any  case,  the  present  process  S2  which  is  commonly 
referred  to  as  ’’stochastic  heterogeneous  Lanchester  square”.  Based 
on  the  structure  of  the  underlying  families  of  assumptions,  this 
designation  seems  to  be  inappropriate. 
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S3a  Heterogeneous  Square  Law  Process 


Assumptions 

1.  There  are  M  Blue  weapon  types  and  N  Red  weapon  types. 

2.  Times  between  shots  fired  by  a  surviving  Blue  type  i  weapon 
are  independent  and  identically  exponentially  distributed  with  mean 
l/r1(i). 

3.  When  a  Blue  type  i  weapon  fires  a  shot  it  is  directed  at  a  Red 
weapon  chosen  from  all  Red  weapons  currently  surviving  according  to  a 
uniform  distribution. 

4.  A  shot  fired  by  a  Blue  type  i  weapon  at  a  Red  type  j  weapon 
is  fatal  with  probability  p1(j,i). 

5.  Assumptions  2-4  hold  for  Red  with  mean  firing  rates  r2(j) 
and  kill  probabilities  p2(i,j). 

6.  The  firing  processes  of  all  weapons  initially  present  are 
mutually  independent. 


Process  Characterization 

Under  the  preceding  family  of  assumptions,  the  survivor  process 

((Bt,  §t))  0  is  a  regular  step  process  with  state  space  fcj  x  _  x  $ 

(M  +  N  time's),  jump  function  X  given  by 


M  N  x. 


y . 


X(x,y)  =  Z  Z  [— P2(i,l)r2(j)y •  +  — ^1  P1(j,i)r1(i)xjL] 
i=l  j=l  J  y 

M  N 

where  x  *  1  =  Z  x,  ,  y  •  1  =  Z  y  ,  transition  kernel  P  given  by 
k=l  k  i=l  1 

M 

y  Z  p1(j,i)r1(i)xi 

P(  (x,y);  (x;  yv...,y^  -  1,  yM))  =  y  .J*1  x(x,y) 

N 

x  E  V2U.,5)r2(j)y. 

P((x,y);  (x1,...,xi  -  Xj^;  y))  =  y-~.\  — - 
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X(x,y) 


and  infinitesimal  generator  Q  given  by 


Q((x,y);  (x; 


•  ? 


yM}) 


y.  M 

y-Ti  Jj  P1(j,i)r1(i)x. 


M  N  x. 

Q((x,y);  (x,y))  =  -  E  E  C— P2(i,  j  )r?(  j )y . 

i=l  j=l  x  x  z  z  3 


+  ~  ;JIP1(j,i)r1(i)xi] 


x.  N 

Q((x,y);  (x^  ..  .jX^-i, .  ..,5^;  y))  =  -  P2(i>j  )r2(j)yj 


For  this  process  the  first-order  approximation  (12)  is  given  by 

x .  N 

E[?tl(?0’  =  (x>  y)]  ~  xi  "  (xTI  2  P2(i,j)r2(j)y.)t 

j=l  J 

provided  x^  >  0,  while  if  y^  >  0 


E[Rj|(B0,  Rq)  =  (x,  y)3  ~  Vj 


y .  M 

(tt4  E  p1(j,i)r  (i)x.  )t; 

y-i  i=1  1  1  1 


here  i  =  1,  . M  and  j  =  1,  . N. 

The  reader  may  well  protest  that  Q  is  of  linear- law  rather  than 
square- law  form,  but  this  really  isn’t  so;  the  factors  (x^)/(x  •  1) 
and  (y . )/(y  •  1)  are  simply  normalizing  constants.  In  view  of  assump¬ 
tion  3,  (y. )/(y  *  1)  is  the  probability  that  a  Blue  weapon,  if  it  fires 
at  a  Red  target  force  y  ,  will  direct  the  shot  at  a  Red  type  j  weapon. 

Implicit  (but  not  really  necessary)  in  the  statement  of  assumption  3 
of  the  Homogeneous  Square  Law  Process  SI  is  the  notion  that  the  one  Red 
weapon  attacked  each  time  a  shot  is  fired  is  uniformly  chosen  from  among 
all  Red  weapons  then  present.  Therefore  the  assumptions  of  Process  S3a , 
and  not  those  of  the  commonly  studied  and  applied  Process  S2,  are  the 
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proper  generalization  of  the  assumptions  of  Process  SI  to  the 
heterogeneous  case.  This  generalization  preserves  dependence  of  mean 
firing  rate  on  only  the  type  of  the  shooting  weapon  and  the  uniform 
fire  allocation.  Hence,  we  believe,  it  is  more  appropriate  to  call 
process  S3 a  a  stochastic  model  of  heterogeneous  Lanchester-square 
attrition  that  it  is  to  so  call  Process  S2. 

The  system  of  differential  equations  analogous  to  this  process  is 
given  by 


✓ 

r  . 
3 


-  r 
— 

E  : 
1=1 


M 

E 

i=l 


k-jCj  j 


i)b . 


b: 


M 

E 

k=l 


N 

E 

j=l 


1^2  ( i  > 


where 

kx(j,  i)  =  p1(j,  i)r1(i) 

and 


k2(i,  j)  =  p2(i,  j  )r2 ( j )  • 


In  particular  this  provides  an  interpretation  of  the  matrices  k  , 
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S3  Heterogeneous  Square  Law  Process 


Assumptions 

The  assumptions  here  are  those  of  the  Process  S3a  but  with 
assumption  3  there  modified  as  follows. 

3.  If  a  Blue  type  i  weapon  fires  a  shot  at  a  time  when  the  Red 
force  is  (the  N- vector)  y  then  the  probability  that  the  one  target 
attacked  is  of  type  j  is  ^(y,]),  j  =  1,  N;  within  the  class  of 

type  j  targets  the  specific  target  attacked  is  chosen  uniformly.  For 
each  i  and  y  ,  ^(y,  •  )  is  a  probability  distribution  over  the 
set  |  1,  . . . ,  N  f  of  Red  weapon  type  indices . 

Similarly,  assumption  5  is  changed  to  include  fire  allocation 
distributions  rij(x,  •  )  for  Red  type  j  weapons  firing  at  a  Blue  force 
of  composition  x  . 

Process  Characterization 

Under  the  family  of  assumptions  given  above,  ((B^,  R^))^^  is  a 
regular  step  process  with  jump  function  X  given  by 


M  N 

\(x,y)  =  E  E  { ^(y,  j)P1(j,i)r1(i)xi 
i=l  j “1 


+  hj(x,i)p2(i,j)r2(j)y.  |  > 


transition  kernel  P  given  by 


P((x,y);  (x;  y±,  ...,  y_.  -  1, 


M 

E^  il»i(y,j)P1(j,i)r1(i)x:L 


N 

E  rij(x,i)p2(i,  j  )r2(j  )y^ 


P((x,  y);  x1, 


•  •  •  3 


and  infinitesimal  generator  Q  given  by 
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Q((x,y);  (x;  yv 


M 

i****^))  =  i|fi(y,j)P1(j,i)r1(i)xi 


Q((x, y);  (x,y))  =  - 


M  N 

S  |  ^(y, j)p1(j,i)r  (i)x.  + 
i=l  j=l 


Tlj(x,i)p2(i,  j  )r2(j  )y .  | 


N 

Q((x,y);  (x1,...,xi  -  1,  ....Xjjjy))  =  ri..(x,i)p2(i,  j  )r2( j )y.j  . 


For  this  process  the  approximation  (12)  has  the  form 


e[bJ| 


N 


Sq)  =  ^x’  y^  ~xi  -Tl  i)P2(i,j)r2(j)y  .It 

Lj  =  j_  J  ^  ^  JJ 


for  i  such  that  x^  >  0.  Similarly,  for  j  such  that  y^  >  0  we  have 

■  r  M  "1 

E[Rjl(B0,  RQ)  =  (x>  ~  yi  ■  E  **(y>  j  )px(j  ,i)r1(i)xiJt  . 


Process  S3a  is  obtained  as  a  special  case  of  Process  S3  by  taking 

y-s 

*i(y’  j)  =  7^1 

and 

x . 

V*’ 

for  all  i  and  j  .  The  generator  Q  of  Process  S3  is  clearly  of 
Lanchester  square  form,  in  the  sense  of  resembling  (5).  Thus,  Process 
S3a  is  a  square  law  process.  This  similarity,  incidentally,  conveys 
additional  information  about  what  the  coefficients  in  (5)  should  mean. 

Effects  such  as  target  priority  and  axiomatically  derived  alloca¬ 
tions  are  included  within  the  large  class  of  S3  processes.  See  the 
Appendix  for  details. 
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We  note  that  the  system  of  differential  equations  analogous  to 
this  stochastic  process  is  given  by 

M 

r'(t)  =  -  E  Y.(r(t),  j)p  (j,  i)r  (i)b.(t) 

J  i=i  1  11 

for  j  =  1,  . ..,  N  and  by 

N 

b:(t)  =  -  E  Y.(b(t),  i)p  (i,  j)r  (j)r  (t) 
j=l  3  1  ^  3 

for  i  =  1,  M.  The  reader  should  be  aware  of  the  potentially 

confusing  notation:  the  r^(i)  and  r^Cj)  are  mean  firing  rates  while 
the  rj(t)  are  surviving  Red  weapons. 

The  assumptions  for  Processes  SI,  S2,  S3  are  all  clearly  similar. 
To  interpret  any  of  these  as  an  area  fire  process  requires  the 
following  reasoning:  since  a  combatant’s  mean  rate  of  fire  is  inde¬ 
pendent  of  the  size  of  the  opposing  force,  he  can  be  envisioned  as 
simply  firing  into  an  enemy- containing  region  according  to  a  Poisson 
process  determined  by  his  own  attributes  (e.g.,  weapon,  location,  ...), 
which  seems  plausible.  But  it  must  then  be  assumed  that  a  shot  kills 
either  exactly  one  enemy  combatant  or  none,  which  doesn't  seem  entire¬ 
ly  reasonable.  Indeed,  we  have  previously  suggested  that  area  fire 
attrition  processes  are  distinguished  by  the  possibility  of  multiple 
kills  arising  from  one  shot.  Hence  it  can  be  argued  that  of  all  the 
processes  presented  here  only  the  Process  A1  is  an  area  fire  process, 
and  that  all  the  others  are  point  fire  processes  which  are  grouped 
into  two  classes  (square-law  and  linear-law)  according  to  the  form 
of  mean  engagement  rates.  For  square-law  processes  the  engagement 
rate  depends  only  on  the  strength  of  the  side  initiating  engagements, 
while  for  linear- law  processes  each  engagement  rate  is  proportional 
to  the  product  of  the  strengths  of  the  two  sides.  These,  of  course, 
correspond  to  rather  different  physical  situations. 
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LI  Homogeneous  Linear  Law  Process 


Assumptions 

1.  All  weapons  on  each  side  are  identical 

2.  The  time  required  for  a  Blue  weapon  to  detect  a  particular 
Red  weapon  is  exponentially  distributed  with  mean  1/s^.  Each  Blue 
weapon  detects  different  Red  weapons  independently. 

3.  A  Blue  weapon  attacks  every  Red  weapon  it  detects;  the 
conditional  probability  of  kill  given  attack  is  q^.  The  attack 
occurs  instantaneously  and  contact  is  immediately  lost.  An  attack 
cannot  occur  without  a  detection. 

4.  Red  weapons  satisfy  the  same  assumptions  with  parameters 

S2’  q2* 

5.  The  detection  and  attack  processes  of  all  weapons  initially 
present  are  mutually  independent. 

Process  Characterization 

Let  k  =  s  q  for  i  =  1,  2. 

Xj  Xj  X) 

Under  assumptions  1  -  5,((Bt,  Rt))  is  a  regular  step  process  with 
state  space  N  x  N,  jump  function  \  given  by 

X(i, j )  =  ij(kx  + 

transition  kernel  P  given  by 

P((i,j),  (i,j  -  D)  = 

P((i,j),  (i  -  l»j))  = 

and  infinitesimal  generator  Q  given  by 


k2) 


k1  +  k2 
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Q((i,j),  (i,j  -  1))  =  k^j 

Q((i,j),  (i» j ))  =  -  ij(^x  +  k2) 

Q((i,j),  (i  -  1,  j)  =  k2ij  . 

As  first  order  approximations  for  this  process  we  have 

ErBt|(B0,  Rq)  =  (i,  j)]  ~  i  -  (kxij)t 

=  i(l  -  kxjt) 

and 

ErRt|(B0,  Rq)  =  (i,  j)]  ~  j  -  (k2ij)t 

=  j(l  -  k2it)  . 

The  similarity  of  Q  and  the  deterministic  homogeneous  Lanchester 
linear  equations  (4)  is  evident,  and  our  process  provides  a  physical 
interpretation  (namely,  inverse  of  mean  time  required  to  kill  a  parti¬ 
cular  opponent)  for  the  coefficients  in  (4). 

The  reader  should  note  the  presence  of  an  assumption  concerning 
detection,  rather  than  firing,  times  in  this  model.  Indeed,  square- 
and  linear-law  models  are  distinguished,  in  our  structure,  in  that 
models  of  the  former  type  include  assumptions  on  firing  times  inde¬ 
pendent  of  the  strength  of  the  opposing  force  while  the  latter  include 
assumptions  on  one-on-one  detection  times. 
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L2  Homogeneous  Linear  Law  Process  With  Engagements 


Assumptions 

1.  All  weapons  on  each  side  are  identical. 

2.  The  time  required  for  a  Blue  weapon  to  detect  a  prescribed  Red 
weapon  is  exponentially  distributed  with  mean  1/s.  .  Each  Blue  weapon 
detects  different  Red  weapons  independently. 

3.  A  Blue  weapon  engages  every  Red  weapon  it  attacks. 

4.  All  engagements  are  binary  (i.e.,  one-on-one).  During  an 
engagement  both  combatants  are  invulnerable  to  other  weapons  and  unable 
to  make  further  detections. 

5.  The  length  of  an  engagement  is  exponentially  distributed  with 
mean  l/u .  An  engagement  terminates  in  exactly  one  of  the  following 
outcomes  with  the  probabilities  indicated,  independent  of  past  history. 


Outcome 

Blue  killed,  but  not  Red 
Red  killed,  but  not  Blue 
Mutual  survival 


Probability 


These  parameters  are  independent  of  the  initiator  of  the  engagement. 

6.  Red  weapons  satisfy  assumptions  2  and  3,  with  mean  detection 
time  1/s 

7.  All  detection  processes,  engagement  lengths,  and  engagement 
outcomes  are  mutually  independent. 


Process  Characterization 

* 

To  study  this  model  we  define  some  further  notation.^  Let  B^_ 
denote  the  number  of  unengaged  Blue  weapons  at  time  t,  R^_  the  number 
of  unengaged  Red  weapons  at  time  t  and  D  the  number  of  binary  engage¬ 
ments  in  progress  at  time  t,  which  is  also  the  number  of  Blue  weapons 
and  Red  weapons  engaged  at  time  t.  We  first  consider  the  three  component 

&  ic 

process  ((Bt,  Rt>  \')')t>0  • 
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Under  assumptions  1  through  7  above,  the  process  ((B*,  R*,  D  ))  q 
is  a  regular  step  process  with  state  space  N  x  N  x  N,  jump  function 
X  given  by 


X(i,j,k)  =  ku  +  ij(s1  +  s2)  , 
transition  kernel  P  given  by 


P((i,j,k),  (i  +  1, j , k 


D) 


kp2u 

ku  +  ij(s^  +  s2) 


P((i,j,k),  (i,  j  +  1,  k 


kp  u 

1))=  V - r~ - 

ku  +  13  (s^ 


) 


P((i,j,k),  (i  +  1,  j  +  1,  k  -  1))  = 


kp3u 


ku  +  ij(s,  + 


s2) 


P((i,j,k),  (i  -  1, 


ij(s1  +  s2) 


^  k  +  1))  ku  +  ij(s1  +  s2) 


and  infinitesimal  generator  Q  given  by 

Q((i,j,k),  (i  +  1,  j,  k  -  1))  =  kp2u 

Q((i,j,k),  (i,  j  +  1,  k  -  1))  =  kpjU 

Q.(  ( i  j  j  j  k ) ,  (i  +  1,  j  +  1,  k  -  1))  =  kp3u 

Q.C  ( i  j  j  5  k ) ,  (i  -  1,  j  -  1,  k  +  1))  =  ij  ( sx  +  s2) 

Q((i,j,k),  (i,j,k))  =  -  [ku  +  ij(s1  +  s2)]  . 

From  the  state  of  i  unengaged  Blue  weapons,  j  unengaged  Red 
weapons,  and  k  engagements  in  progress,  four  states  may  be  reached, 
corresponding  to  the  events  ’’termination  of  one  engagement  with  kill 
of  a  Red  weapon”,  ’’termination  of  one  engagement  with  kill  of  a  Blue 
weapon”,  ” termination  of  one  engagement  with  mutual  survival”,  and 
initiation  of  one  additional  engagement”.  These  occur  at  the  rates 
indicated  by  the  expressions  for  Q  (which,  as  usual,  assume  that 
i  >  0,  j  >  0,  k  >  0  and  must  be  suitably  adjusted  for  other  states). 

Based  on  its  assumptions  and  the  term  in  Q  expressing  engage¬ 
ment  initiation,  this  process  is  of  linear  law  type.  But  our 
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principal  interest  is  in  the  survivor  process  ((Bt,  Rt))  defined  by 


Bt  =  B?  + 


and 


R 


t 


+  D. 


since  the  number  of  Blue  weapons  surviving  at  time  t  is  the  sum  of  the 
numbers  of  unengaged  and  engaged  Blue  weapons  surviving  at  time  t,  and 
similarly  for  Red.  The  process  ((Bt,  R  ))  is  not,  unfortunately,  a 
regular  step  process  and  we  have  not  yet  characterized  it  in  a  useful 
manner.  However,  expectations  of  B^  and  R  are  easy  to  compute  from 
those  of  B*  j.  R*>  and  D^.,  because 

E[Bt]  =  E[B*]  +  E[Dt] 

and  E[R  ]  =  E[R*]  +  ErDt]  . 

The  latter  expectations  might  be  computed,  since  ((B*,  R*,  D  ))  _ 

L  L  L  L  ^  U 

is  a  regular  step  process,  using  the  methods  described  in  Section  IV. 

To  compute  first-order  approximations  to  E[Bt]  and  E[Rt]  we  first 
note  that  ( 12)  and  the  form  of  the  generator  Q  of  this  process  imply 
that 

E[B*| (Bg,  RJ,  DQ)  =  (i,j,k)]  ~  i  -  [ij(s1  +  s2)  +  ku(p2  +  p3)]t  , 

that 

E[R*|  (Bg,  R*,  Dq)  =  (-  ,  j ,k)]  ~  j  -  rij(S;L  +  s2)  -  ku(P;L  +  p3)]t  , 

and  also  that 

E[Dt|(B*,  R*,  D0)  =  (i,  j  ,k)3  ~  k  -  [ku  -  ij(s;L  +  s2)]t  . 
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Thus 


ECBtl  (B0’  R0’  D0}  =  (i»J,k)] 

=  E[B*|(B*,  R*,  Dq)  =  (i,j,k)]  +  E[Dt|(B*0,  R*,  DQ)  =  (i,j,k)] 
~  i  -  ij(s1  +  s2)t  +  ku(p2  +  p3)t  +  k  -  kut  +  ij(s1  +  s,,)t 
=  i  +  k  -  (p  ku)t  ; 

we  similarly  obtain 

E[Rt|(Bj,  Rq,  Dq)  =  ( i , j , k ) ]  =  j  +  k  -  (p2ku)t  . 

The  assumption  that  the  two  weapons  in  a  binary  engagement  be 
invulnverable  to  other  opposition  weapons  is  questionable  (an  engaged 
weapon  would  seem  to  be  relatively  more  vulnerable  than  an  unengaged 
weapon)  and  alternatives  should  be  sought. 
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L3  Heterogeneous  Linear  Law  Process 


Assumptions 

1.  There  are  M  Blue  weapon  types  and  N  Red  weapon  types. 

2.  The  time  required  for  a  Blue  type  i  weapon  to  detect  a  parti¬ 
cular  Red  type  j  weapon  is  exponentially  distributed  with  mean 
l/s^jji).  Each  Blue  weapon  detects  different  Red  weapons  independently. 

3.  A  Blue  weapon  attacks  every  Red  weapon  it  detects.  The  condi¬ 
tional  probability  that  a  Blue  type  i  weapon  kills  a  Red  type  j  weapon, 
given  detection  and  attack,  is  p^(j,i). 

4.  Red  weapons  satisfy  assumptions  2  and  3  with  mean  detection 
times  l/s2(i,j)  and  kill  probabilities  p2(i,j)  describing  Red  type  j 
weapons  opposing  Blue  type  i  weapons. 

5.  Detection  and  attack  processes  of  all  weapons  present  are 
mutually  independent. 

Process  Characterization 

For  q  =  1,  2  and  appropriate  k  and  £  define  k^(k,£)  = 

Sq(k, £)Pq(k, £) .  Then  under  assumptions  1  through  5  above,  the  vector- 
valued  survivor  process  ((B  ,  R  ))  Q  is  a  regular  step  process  with 
state  space  x  ...  x  N  (M  +  N  times’),  jump  function  X  given  by 


M  N 

\(x,y)  =  E  E  x.y . [k.( j , i)  +  k2(i,j)]  , 

i=l  j=l  J 


transition  kernel  P  given  by 


M 

y-j  E  vi(j’i)xi 

J  *5=1 


P((x,y);  ( x 5  y^>***>yj 


N 

x±  I  k2(i,j)y^ 


P(  ( Xj  y )  >  ( x^>  •  • .  >  x_^  1 
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and  infinitesimal  generator  Q  given  by 


Q((x,y);  (x;  y±, ... ,y ^ 


,yN»  =  Vj 


M 

E  k,(j,i)x. 
i=l  J 


M  N 

Q((x,y);  (x,y))  =  -  E  E  x .y.[k ,(j,i)  +  k  (i,j)] 

i=l  j=l  3  z 


N 


Q( ( x, y ) }  ( xx, • • . 5 xi  1, . . . , x^ j y ) )  x^  E  k2(i5j)y-  . 

j=l  J 

Expectations  of  numbers  of  survivors  in  this  process  may  be 
approximated  using  (12)  in  the  following  manner: 


ECBtlCB0,  R0)  =  (x,y)] 


N 


and 


E[R^|(B0,  Rq)  =  (x,y)] 


x±  -  [x.  E  k2(i,j)yj]t 
j=l 

N 

x  [1  -  t  E  k  (i,j)y  )] 

i=l  3 


M 

y  -  [y.  E  k  (j,i)x  It 
3  3  i=l  1 


M 


y  .[1  -  t  E  k,  ( j  ,i)x. )]  , 

3  i=l  ±  1 


for  i  =  1,  . M  and  j  =  1,  N. 


Based  on  its  underlying  assumptions  this  process  is  clearly  the 
heterogeneous  analog  of  the  stochastic  homogeneous  linear  law  Process 
LI.  Its  infinitesimal  generator  Q  is  similar  to  the  system  (6), 
which  comfirms  that  (6)  is  an  appropriate  deterministic  heterogeneous 
linear  model  and  provides  the  proper  interpretations  of  the  coeffi¬ 
cients  in  (6);  namely,  the  inverse  of  the  mean  time  required  to  detect 
and  kill  a  particular  opposition  weapon  of  a  given  type. 
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Ml  Heterogeneous  Mixed  Law  Process 


Assumptions 

1.  Each  side  possesses  exactly  two  types  of  weapons. 

2.  Times  between  shots  fired  by  a  surviving  type  1  Blue  weapon 
are  independent  and  identically  exponentially  distributed  with  mean 
l/vr 

3.  Each  such  shot  is  directed  at  a  target  chosen  according  to  a 
uniform  distribution  on  the  set  of  all  Red  weapons  surviving  at  that 
time.  If  the  target  is  a  Red  type  j  weapon  it  is  killed  with  proba¬ 
bility  p-^j). 

4.  Red  type  1  weapons  satisfy  assumptions  2  and  3  with  mean 
firing  rate  r^  and  kill  probability  p2(i)  against  a  Blue  target  of 
type  i. 

5.  The  time  required  for  a  Blue  type  2  weapon  to  detect  a  parti¬ 
cular  Red  weapon  of  type  j  is  exponentially  distributed  with  mean 
l/s1C j ) .  A  fixed  Blue  type  2  weapon  detects  different  Red  weapons 
independently  of  one  another. 

6.  A  Blue  type  2  weapon  attacks  every  Red  weapon  it  detects. 

The  conditional  probability,  given  detection  and  attack,  that  a  Red 
weapon  of  type  j  is  killed,  is  q-^j).  The  attack  occurs  instan¬ 
taneously  and  contact  is  lost. 

7.  Red  type  2  weapons  satisfy  assumptions  5  and  6  with  mean 
detection  times  l/s2(i)  and  kill  probabilities  q2(i)  against  type  i 
weapons . 

8.  All  weapons  detect  and  attack  independently  of  one  another. 


Process  Characterization 


Under  assumptions  1  through  8  above,  the  survivor  process 
((B1.  B~ 


Jt’ 


Rt> 


2 

R  ))  is  a  regular  step  process  with  state  space 
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N  x  N  x  N  x  N,  jump  function  \  given  by 


X(i, j, 

transition  kernel 

P((i> 3  >k,jO; 

P((i, j ,k,£); 

P((i, j,k,£); 


p  (l)k  +  p A  2)1 

c,£)  =  ir,  ( - r - - - ) 

lv  k  +  £  J 

+  3  [ks1(l)q1(l)  +  £s1(2)q1(2)] 

Po(l)i  +  Po(2)j 

+  kr0  (-i-, - - — - - ) 

2  v  1  +  3  J 

+  £  [is2(l)q2(l)  +  js2(2)q2(2)] 
3  given  by 


ir  P,(l)  Y~n  +  3^  (l)qAl) 


;i,j,k,£  -  1))  = 


:i  -  l,j,k,£))  = 


;i,j  -  l,k, £))  = 


iripi(2)  FTT  +  3^siC2)q1(2) 
Ui,j,k,£) 

kr2p2(l)  j  ±  ,  +  £is2(l)q2(l) 
X(  fjj  £ ) 

kr2p2(2)  j  j-  .  +  £js2(2)q2(2) 
\(i, j,k,£) 


and  infinitesimal  generator  Q  given  by 


Q((i,j,k,£);  (i,j,k  -  1,£))  =  ir^Cl)  -  +  jks  (i)q  (i) 

Q((i,  j,k,£);  (i,j,k,£  -  1))  =  ir1P1(2)  -R  ~  -  +  j£s1(2)q,(2) 

p,(l)k  +  p,(2)£ 

(i,j,k,£))  =  -  |  ir^— - - ) 

+  j[ks1(l)qi(l)  +  £s1(2)q1(2)] 


p2(i)i  +  PpC  2 ) j 
+  kr0  (-£ - ^—4 - ) 


a-  +  3 


+  £[is2(l)q2(l)  +  js2(2)q2(2)]} 


Q((i,j,k,£);  (i  -  l,j,k,£))  =  kr2p2(l)  + 


Q((i,j,k,£);  (i,j  -  l,k,£))  =  kr2p2(2)  j-j-j  + 


For  the  approximation  (12)  we  have 
E[bJ|(bJ,  Bq,  rJ,  Rq)  =  (i,j,k,i)] 


£is2(l)q2(l) 
£js2(2)q2(2)  . 


~  i  -  tkr2P2(l)  ~  £is2(l)qs(l)]t 
kr  p  (1) 

=  i(l  -  tr-T£_£^_  is2(i)q2(l)])  , 

E[B^| (Bq,  Bq,  Rq,  Rq)  =  <i,j,k,£)] 
kr  p  (2) 

~  j(l  -  t[  jLZj  +  £s2(2)q2(2)])  , 
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E[RJ|(B0>  Bq,  Rj,  Rq)  =  (i,j,k,i)] 
ir  p  (1) 

~  R(1  -  tr-  ■  y  ±  J  +  js^Dq^l)])  , 

and 

e[R^|(bJ,  b*, 

~  Jt(l  - 

respectively. 

The  four  states  which  can  be  reached  from  state  (i,j,k,jO  correspond 
to  deaths  of  the  four  different  weapon  types  and  are  entered  at  the 
’’rates”  indicated  in  the  expression  for  Q  .  Comparing  Q  here  with 
the  generators  of  the  Processes  SI  and  L3  yields  the  ”mixed  law” 
interpretation  of  this  process,  as  does  a  comparison  of  the  families 
of  assumptions  underlying  the  three  processes. 

No  heterogeneous  mixed  laws  have  been  discovered  before.  Without 
families  of  assumptions  to  use  in  making  generalizations,  it  was  not 
known  what  appropriate  process  would  be.  Moreover,  the  Process  Ml 
can  be  extended  to  allow  an  arbitrary  number  of  weapon  types  on 
each  side. 


R0’  V 

xl'p1(2) 

tc-rVr-"  jsi(2)qi(2)^  > 
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Mia  Homogeneous  Mixed  Law  Process 


Assumptions 

1.  All  weapons  on  each  side  are  identical. 

2.  Red  weapons  satisfy  assumptions  2  and  3  of  process  Ml  with 
parameters  r,  p. 

3.  Blue  weapons  satisfy  assumptions  5  and  6  of  process  Ml  with 
parameters  s,  q. 

4.  Detection  and  attack  processes  of  all  weapons  are  mutually 
independent. 

Process  Characterization 

Under  assumptions  1  to  4,  the  process  ((B^,  R^)  )1_>g  is  a  regular 
step  process  with  state  space  N  x  N,  jump  function  \  ""given  by 

\(i,j)  =  irp  +  ijsq  , 

transition  kernel  P  given  by 


P((i,j),  (i,j  "  D)  = 


rp 


rp  +  jsq 


P( (i ? j  ) ?  (i  -  l,j))  =  rp3vq-j-s- 

and  infinitesimal  generator  Q  given  by 

Q((i,j),  (i,  j  -  D)  =  irp 

Q((i,j),  ( ij j ) )  =  "  (irP  +  ijsq) 

Q(  ( i  ?  j  )  ?  (i  "  1,  j))  =  ijsq  . 

First-order  approximations  in  this  case  are  given  by 
E[Bt|(B0,  Rq)  =  (i, j )]  =  i  "  jrpt 

and 

E[Rt|(B0,  Rq)  =  (i,j)3  =  jC1  -  t[isq]), 

respectively. 
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APPENDIX 


PROBABILISTIC  TECHNICALITIES  AND  PROOFS  OF  RESULTS 

We  give  here  the  proofs  of  the  characterization  theorems  in  the 
main  body  of  the  paper.  First,  however,  some  probabilitistic  concepts 
related  to  regular  step  processes  will  be  discussed  in  further  detail. 

Let  E  be  a  countable  set  with  discrete  c-algebra  E  .  A  Markov 
kernel  on  the  measurable  space  (E,  E)  is  a  mapping  P  of  E  x  E  into 
[0,  1]  such  that  A  -*  P(i,  A)  is  a  probability  measure  on  £  for  each 
i  e  E.  Since  E  is  discrete  the  probability  P(i,  •)  is  determined 
by  its  values  on  singleton  sets  and  P  may  hence  be  considered  as  the 
matrix  P  defined  by 

P(i,  j )  =  P(i,  |  j} )  . 

If  f  is  a  bounded  or  nonnegative  function  on  E  ,  Pf  denotes  the 
function  defined  by 


Pf (i)  =  E  P(i,  j)f(j)  • 
jeE 

Given  a  function  X:  E  -  [0,  ®)  and  a  Markov  kernel  P  on  (E,  |) 
such  that  P(i,  i)  =  0  for  all  i  e  E  such  that  X(i)  >  0,  there  exists 
a  continuous  parameter  Markov  process  X  =  (Q,  £,  Xt,  P1)  with  state 
space  (E,  E)  satisfying  the  following  intuitive  description.  If  the 
process  enters  a  state  i  with  \(i)  =  0  it  remains  there  forever 
after  (such  states  are  said  to  be  absorbing).  When  the  process  enters 
a  state  i  with  X(i)  >  0  its  sojourn  time  there  is  exponentially 
distributed  with  mean  1/X(i)  and  independent  of  the  past  history  of 
the  process.  At  the  end  of  this  time,  the  process  jumps  (instantan¬ 
eously)  to  a  new  state  of  E  according  to  the  probability  distribution 
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P(i,  • ),  independent  of  the  present  time,  the  length  of  its  sojourn 
in  state  i  and  all  previous  jumps.  Hence,  successive  states  entered 
form  a  Markov  chain  with  transition  matrix  P  . 

DEFINITION.  X  is  the  regular  step  process  with  jump  function  X 
and  transition  kernel  P  . 

P1  denotes  the  probability  law  of  the  process  given  that  Xq  =  i; 
expectation  with  respect  to  this  measure  is  denoted  by  E1,  so 

E1[Y]  =  J  Y(uu)P1(duj) 

n 

for  suitable  F-measurable  random  variables  Y. 

Let  bE  denote  the  set  of  all  bounded  (and  necessarily  E-measurable) 
functions  on  E  .  Then  the  relations 

Ptf(i)  =  E1[f(Xt)] 

define  a  semigroup  (P^)  q  of  bounded  linear  operators  on  bE,  called 
the  transition  function  of  X  .  If  f  is  the  indicator  function 
l|j|  of  the  singleton  set  -{jj  then 

Ptf(i)  =  P1{Xt  =  j (  s  Pt(i,  j) 


and  for  any  g  , 


Ptg(i)  =  s  Pt(i,  j)g(j)  . 
j 

Hence  (P  )  is  represented  by  a  family  of  stochastic  matrices;  we  do 
not  notationally  distinguish  the  operator  P^  and  the  corresponding 
matrix.  Pq  is  the  identity  matrix. 

The  transition  function  (P  )  is  uniquely  determined  by  its 
infinitesimal  generator,  which  is  the  (possibly  unbounded)  linear 
operator  Q.  defined  by 


56 


Qf  (i)  =  lim 
1 10 


Ptf(i)  -  f(i) 
t 


Q  is  represented  by  the  matrix  Q  given  by 

Pt(i>  j)  -  6(i,  j) 

Q(i>  j)  =  lim  - r -  , 

ti  0 

since  Pj-  is  the  identity  matrix. 

Stated  alternatively,  the  matrix-valued  function  t  -♦  P  is 
differentiable  (componentwise)  and 


for  all  t  ;  in  particular 


The  transition  semigroup  is  given  in  terms  of  the  infinitesimal 
generator  Q  by  formal  solution  of  the  system  P'  =  QP  of  linear 
differential  equations  by  exponentials;  thus 


P 


t 


e 


tQ 


where 


.«  =  l  ^ 

n=0  * 


For  regular  step  process,  Q  is  given  explicitly  in  terms  of 
the  jump  function  \  and  transition  kernel  P  by  the  relations 


Q(i,  J)  = 


-  X(i) 
\(i)P(i,  j) 


if  j  =  i 
if  j  7*  i  > 


which  may  be  proved  computationally. 
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Fundamental  to  our  claim  that  the  stochastic  processes  characterized 
in  the  following  theorems  are  analogs  of  Lanchester ’ s  differential 
equation  models  of  combat  is  the  ’’differential”  interpretation  of  the 
infinitesimal  generator  Q  .  From  the  representation  P  =  e^  we  see 
that 


Ph(i,  j)  =  6(i,  j)  +  hQ(i,  j)  +  o(h) 

as  h  -  0.  In  particular  for  j  /  i 

Ph(i,  j)  =  hQ(i,  j)  +  o(h)  ,  h  -  0  . 

That  is,  the  probability  of  a  jump  from  i  to  j  in  the  interval 
(0,  h]  is  approximately  h  •  Q(i,  j)  so  that  one  can  interpret 
Q(i,  j )  as  the  infinitesimal  rate  at  which  the  process  moves  from 
state  i  to  state  j  . 

We  next  present  proofs  of  the  characterization  theorems  for  each 
of  our  stochastic  attrition  processes. 
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A1  Homogeneous  Square  Law  Area  Fire  Process 

We  construct  this  process  in  the  following  manner.  Let  E  =  N  x  U, 
where  N  =  {  0,  1,  2,  . ..},  the  set  of  natural  numbers.  As  the  sample 
space  Q  we  take  the  family  of  all  functions  id  =  (id^,  ujg)  mapping 
[0,  oo )  into  E  with  the  properties  that 

a)  u)  is  right  continuous  on  [0,  00 );  for  each  time  t  >  0, 

lim  u)(u)  =  u)(t); 
uit 

all  limits  in  E  are  taken  in  the  product  of  the  discrete  topologies 
on  the  factor  spaces; 

b)  id  possesses  left-hand  limits  on  (0,  oo);  that  is,  for  each 
t  >  0 


lim  f(u) 
u  T  t 

exists . 

We  define  two  families  (B.  and  (R.  ).  n  of  coordinate  random 

t  t>0  t  t>0 

variables  on  n  by 

Bt(  id)  =  u^Ct)  ,  t  >  0 

and 

Rt(uj)  =  u>2(t)  ,  t  >  0  , 

respectively.  The  interpretation  is  that  Bt  is  the  number  of  sur¬ 
viving  Blue  combatants  at  time  t  and  R  the  number  of  Red  combatants 
surviving  at  time  t. 

Further,  we  let  £  be  the  history  generated  by  the  random 

variables,  {(B  ,  R  )  :  0<s<t},  which  can  be  interpreted  as  the 
s  s 

history  of  the  attrition  process  up  until  time  t,  and  let 
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F  =  a((B  ,  R  )  :  u  >  0), 

_  ^ '  u’  u'  —  ' 

which  is  the  entire  history  of  the  process. 

For  each  (i,  j)  e  E,  let  denote  the  probability  law  of  the 

attrition  process  governed  by  assumptions  1-5  on  page  23,  conditioned 
on  the  event  {(Bq,  Rq)  =  (i,  j)}.  While  one  can  give  the  details  of 
the  construction  of  these  measures  (beginning  with  probabilities 
assigned  to  appropriate  cylinder  sets  and  proceeding  through  an  appli¬ 
cation  of  the  Kolmogorov  Extension  Theoren)  such  an  approach  is  not 
appropriate  in  this  exposition.  So,  instead,  we  take  for  granted  the 
existence  of  such  probabilities  and  derive  in  our  Theorems  characteri¬ 
zations  of  the  resultant  stochastic  attrition  process. 

(1)  THEOREM.  Subject  to  assumptions  1-5  of  the  family  Al,  the 
process 

(0,  E,  Et>  (Bt,  Rt),  P(l,:,)) 


is  a  regular  step  process  with 

a)  state  space  (E,  E);  we  remind  that  E  =  N  x  N  and  E  is  the 
discrete  a-algebra; 

b)  jump  function  \  given  by 


(2)  \(i,j)  = 


ir^l  -  (1-P1)j]  +  jr2[l  -  ( 1“P2 )1]  if  i  >  0,  j  >  0 
0  if  i  =  0  or  j  =  0 


t)  transition  kernel  P  given  for  states  (i,j)  with  i  >  0  and 
j  >  0  by 
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(3) 


P((i,j);  (i,4)) 


p 


1-3 

1 


y 


o  <  i  <  j 


(i  -  p2) 
X(i, j ) 


k 


0  <  k  <  i  ; 


d)  infinitesimal  generator  Q  given  for  states  (i,j)  with  i  >  0 
and  j  >  0  by 


Q((i,j);  (i, £))  =  -  P±)1  ,  0  <  £  <  j 

Q((i,j);  (i, j ))  =  -  (ir^Cl  -  ( l-p1)^i 3  +  jr2[l  -  (l-Pj)1]) 

Q((i,j);  (k,j))  =  jr2(y(l  -  Px)k  p^"k  ,  0  <  k  <  i 


If  i  =  0  or  j  =  0,Q((i,j);  x)  =  0  for  all  x  e  E. 

PROOF.  For  detailed  arguments  we  refer  to  the  proof  of  Theorem  (6). 

First  of  all,  note  that  if  some  Blue  combatant  fires  a  shot  at  an 
instant  when  there  are  j  surviving  Reds,  the  probability  that  l 
Reds  survive  is  the  binomial  probability 

e>‘  -  -r  • 

in  particular,  the  probability  that  the  shot  causes  one  or  more 
fatalities  is 

q-L  =  i  -  (i  -  p1)^  • 


Thus  if  we  denote  by  Sn  and  S_  the  times  of  the  first  fatality- 
causing  shots  fired  by  Blue  and  Red,  respectively,  then  with  respect 
to  the  probability  Sg  and  SR  are,  each  conditioned  on  non¬ 

occurrence  of  the  other,  independent  and  exponentially  distributed 
with  means  1/ir^q^  and  l/jr2q2,  respectively.  It  follows  (see  the 
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proof  of  Theorem  (6))  that  the  time 


V  =  infjt:  (Bt,  Rt)  *  (BQ,  RQ)\ 


of  the  first  change  of  state  of  the  process  is  exponentially  distributed 
with  mean  l/(ir1q1  +  j r'2Cl2 ^  (with  resPect  to  of  course;  note  also 

that  is  a  function  of  j  and  q^  a  function  of  i  ,  even  though  we 
have  suppressed  this  dependence  in  our  notation)  and  that,  moreover 


p<1,])|<BV+>  V5  e  l(i’  l): 
=  P(iJ)|sB  <  SR| 


“A 


ir'l<^l  +  ~*r2<^2 


0  <  l  <  j  }f 


while 


P^’j)|(Bv+>  Ry+)  €  |(k,  j):  0  <  k  <  ill 


=  P(i’j)|SR  <  sBl 


jr2q2 


lrlql  +  jr2q2 


Let  K-  be  the  number  of  fatalities  caused  by  the  first  fatality- 
B 

causing  shot,  if  there  is  one,  fired  by  a  Blue  combatant.  Then 
conditioned  on  the  event  I SB  <  SR|,  Kg  is  "binomially"  distributed  on 

{ C,  ...,  j  -  lj,  with  respect  to  P^’^,  .in  the  sense  that 


?(!> j )|v  =  £js 


a  =  o5 


j  -  i 
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the  normalization  is  required  in  order  that 


j-1 


E  P^’^lKg  =  *|SR  <  Sp|  =  1  . 


£=0 


B  R' 


We  then  have  for  0  <  SL  <  j 

‘>(i,3)KBv+>  Rv+)  =  (i>  1)1 
=  f(i,;i)lSB  <  v  KB  =  j  - 

=  P(1>;i)tSB  <  SR)P<i>j)|KB  =  j  -  £|Sb  <  SRt 

lrA  PB'f 

irA  +  3r2q2  1  -  <l-p 


ir  Q 
141 


iriQl  +  :r2q2 


 iri(i)(i- 


p/  p3-‘ 


*X(i,  j) 


=  P((i,j),  (i,£)) 
where  P  is  defined  by  (3)  . 

We  leave  to  the  reader,  should  he  desire  the  details,  the  entirely 
analogous  proof  that 

P(i’jH(Bv+,  Ry+)  =  (k,j)| 

i  , .k  i -k 
Jr2  (  X-Pg)  P2 

=  P((i,j)j  (k, j ) ) 


for  0  <  k  <  i  . 
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Hence  with  respect  to  P^’  ,  the  random  variable  V  is 
exponentially  distributed  with  mean  l/\(i,j)  where  \  is  defined  by 
(2)  and  (By+,  R^+)  is  distributed  as  P((i,j);  •).  The  extension  to 
arbitrary  fixed  times  is  omitted.  n 

The  binomial  distributions  postulated  in  assumptions  3  and  4  have 
the  (perhaps)  unrealistic  effect  of  permitting  each  side  with  positive 
probability  to  annihilate  the  other  with  a  single  shot. 

It  is,  from  a  theoretical  standpoint,  quite  acceptable  to  replace 
the  family  of  binomial  distributions  by  two  families  JcPjC*):  j  >  1} 
andjY.(-):  i  >  l}  of  probabilities  on  the  nonnegative  integers  with 
the  property  that  for  each  j  ,  cd_.  is  a  probability  distribution  on 
jo,  ...,  j[  and  for  each  i  ,  y.  is  a  probability  distribution  on 
jo,  . i]  .  The  interpretation  is  then  that  cp^  ( X )  is  the  proba¬ 
bility  that  a  shot  fired  by  a  Blue  combatant  at  an  instant  when  there 
are  j  surviving  Reds,  kills  exactly  j  -  SL  of  those  j  Reds,  and 
similarly  for  the  Y. .  We  may  then,  with  only  notational  changes,  obtain 
the  following  generalization  of  Theorem  (1). 

(  4)  THEOREM.  The  process 

(0,  E,  Et,  (Bt,  Rt),  P^’35) 


is  a  regular  step  process  with 

a)  state  space  (E,  E); 

b)  jump  function  \  given  by 


X(i,  j)  = 


irx(l  -  cpj(j))  +  Z)>2(1  -  Y±(i) ) 


if  i  >  0,  j  >  0 
if  i  =  0  or  j  =  0; 


/  0 

l 
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c)  transition  kernel  P  given  for  states  (i,j)  with  i  >  0  and 
j  >  0  by 


P( 


ir  cp-(A) 

o<«<j 

jr2  ^(k) 


P((i,j);  (k,j))  = 


X(i,j ) 


0  <  k  <  i  ; 


d)  infinitesimal  generator  Q  given  for  states  (i,j)  with  i  >  0, 
j  >  0  by 

Q((i,j),  (i,X))  =  irx  cp^- ( A )  ,  0  <  l  <  j 

Q(  (  i  j  j  )  >  (ij  j  )  )  =-rir1(l  -  cPj(j))  +  jr2(l  -  Y±(  i)  )3 

Q((i»j),  (k,j))  =  jr2  Y±(k)  ,  0  <  k  <  i  . 

If  i  =  0  or  j  =  0,  Q((i,j),x)  =  0  for  all  x  e  E. 


There  is  no  need  to  repeat  the  proof  of  (1);  it  goes  through 
here  with  only  notational  changes. 

One  pays  a  price,  however,  for  such  generality  and  that  price  is 
that  most  such  distributions  jcp^l  and  do  not  arise  from  a  set  of 

physical  assumptions  concerning  the  attrition  process  being  modeled. 
The  binomial  distributions,  for  all  their  other  unpalatable  character¬ 
istics,  at  least  do  have  a  physical  basis  in  the  assumption  that  each 
shot  fired,  for  example,  at  a  Red  target  force  of  j  survivors  kills 
each  Red  with  probability  p^^  and  different  Reds  independently  of  one 
another,  for  then  the  number  of  Reds  surviving  the  shot  is  binomially 
distributed  with  parameters  (j,  1  -  p^). 

If  the  distributions  |cpjf  and  |y^(  are  all  of  the  form 
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CPj(j-l)  =  P-j_ 

Cpj(j)  =  1-P]_ 

(5) 

Yi(i-1)  =  P2 
Tid)  =  1-P2  , 

where  0  <  <  1,  0  <  p0  <  1,  then  the  Process  SI  below  is  obtained. 
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SI  Homogeneous  Square  Law  Process 


The  canonical  stochastic  process 


(0,  E,  Et,  (Bt,  Rt)) 

remains  unchanged.  The  relevant  assumptions  are  the  family  of 
assumptions  1  through  5  on  page  27.  At  the  risk  of  possible  confusion, 
but  in  order  to  prevent  a  hopeless  proliferation  of  notations  for 
various  families  of  probability  measures  on  (Q,  F)  we  again  denote  by 
p(i,j)  the  probability  law  of  the  attrition  process  governed  by  this 
family  of  assumptions,  subject  to  the  initial  conditions  Bg  =  i,  Rg  =  j 
Thus,  while  the  family  of  assumptions  has  changed,  the  notation  has  not 
the  reader  should  keep  this  in  mind  here  and  throughout  the  exposition. 


(6)  THEOREM.  Under  assumptions  1-5  of  the  family  SI,  the  process 

CO,  E,  Et,  (Bt,  Rt),  P(l,3)) 


is  a  regular  step  Markov  process  with 

a)  state  space  (E,  E),  where  E.  is  the  a-algebra  of  all  subsets 
of  E; 

b)  jump  function  X  given  by 


(7) 


X(i,  l)  = 


if  i  >  0  and  j  >  0 
if  i  =  0  or  j  =0 


c)  transition  kernel  P  given  for  states  (i,  j)  e  E  with  i  >  0 
and  j  >  0  by 


(8) 


P((i>  j  “  1))  ic^  +  jc^ 

jc2 

j)>  (i  "  3  ) )-  To"  a.  -ic^ 
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d)  infinitesimal  generator  Q  given  for  states  (i,  j)  with 
i  >  0  and  j  >  0  by 

Q((i,  j),  (i,  j  -  1))  =  ic-L 

(9)  Q((i,  j),  (i,  j))  =  -  (icx  +  jc2) 

Q((i>  l),  (i  -  1,  j))  =  jc2  . 

If  i  =  0  or  j  =0,  Q((i,  j),  x)  =  0  for  all  x  e  E. 

PROOF.  Let  Eq  =  |(i,  j)  :  i  >  0,  j  >  0[.  Since  there  clearly  is  no 

further  attrition  once  one  side  is  annihilated,  every  state  in  E  -  E 
is  absorbing.  Also,  with  probability  one  the  state  (0,  0)  is  never 
entered  unless  the  process  begins  there. 

We  begin  by  showing  that  if  (i,  j)  e  Eg  and  t  >  0,  then 

p(iJ){(Bu>  ru)  =  (i,  j)  for  all  u  e  TO,  tl} 

=  P(l’D)|(Bt,  Rt)  =  (i,  j)} 

=  exp  r  -  (ic1  +  jc2)tl  . 

Let  V  denote  the  time  at  which  the  process  first  changes  state, 
which  in  this  case  is  also  the  first  time  at  which  a  fatal  shot  is 
fired.  We  are  trying  to  compute  p^’-^jy  >  t} . 

Fix  attention  on  a  particular  Blue  combatant,  say  the  kth  one  and 
consider  the  times  H  o,  bS2^’  at  which  he  w°tld 

fire  shots  if  he  survived  and  continued  firing  forever.  By  assumptions 
2)  and  3)  the  sequence  CBs^^)n>Q  is  a  simPle  Poisson  process  with 
rate  r^  .  Suppose  in  addition  "that  all  shots  fired  by  all  combatants 
have  only  "hypothetical”  effects  in  the  sense  that  no  combatants  are 
killed  but  a  record  is  kept  of  shots  which  would  have  been  fatal  with¬ 
out  this  restriction.  This  converts  the  attrition  process  into  a  non¬ 
terminating  process  in  which  each  combatant  continues  firing  shots 
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forever,  some  of  which  cause  "hypothetical"  kills  of  opponents.  The 
hypothetical  process  not  only  has  nicer  behavior  than  the 
original  attrition  process  but  also  serves  to  describe  the  attri¬ 
tion  process  because  of  the  following  two  properties: 

(a)  The  time  of  the  first  hypothetical  kill  is  the  time  at  which 
the  attrition  process  first  changes  state; 

(b)  The  side  suffering  the  first  hypothetical  casualty  also  suffers 
the  first  casualty  in  the  attrition  process. 

Thus  consideration  of  the  process  of  hypothetical  kills  suffices 
by  (a)  to  compute  the  jump  function  of  the  attrition  process  and  by 
(b)  to  compute  its  transition  kernel,  as  well  as  by  (a)  and  (b) 
together  to  establish  the  Markov  property.  Computations  involving 
the  process  of  hypothetical  kills  are  easy  to  do,  so  we  shall  use 
throughout  the  approach  of  considering  processes  of  hypothetical  kills. 

By  the  Random  Sampling  Theorem  for  Poisson  processes,  the 
times  (BS^)n>0  of  hypothetically  fatal  shots  fired  by  the  k1"^  Blue 
combatant  form- a  Poisson  process  with  rate  c^  =  r^p^. 

Similarly,  for  each  i  ,  the  set  (  S^^)  of  times  of  hypothetically 

t  fl  ^  ^ 

fatal  shots  fired  by  the  £,  Red  combatant  is  a  simple  Poisson  process 
with  rate  c^. 

By  assumption  5)  the  random  variables  ...,  nS~^, 

d  L  B  1 

^S^1^  ,  ....  -.S^^  are  mutually  independent  with  respect  to  P^’-^. 

K  1  R  -L 


Furthermore,  V  >  t  if  and  only  if  ^S^  J  >  t  for  k  =  1 


•  •  •  9 


i  and 


qU) 

R^  >  t  for  SL  =  1,  j,  since 


V  =  m: 


B°1  >  R°1  >  *••»  R’ 


(i)  e(l) 


,S 
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Therefore 


pC1)  D  )  j  v  >  t}=  P(l’:,)|BS^k)  >  t,  k  =  1,  ...,  i;  rs[Z>)  >  t,  JL  =  1,  .  ..,  j 

=  p^«Jbs?>  >  t}J  ^  n  p(1,;5J|rs^j  >  t}J 


I"  n 

Lk=l 


-C.  t-1 


j  -c0t 


n  e 
JU=l 


=  exp  r  -  (ic1  +  jc2)t]  . 

The  same  reasoning  and  %emory  less"  properties  of  the  exponential 
distribution  apply  to  show  that  for  each  t  and  s 

P(1’3){(Bu,  Ru)  =  (Bs,  Rs)  for  all  u  e  [s,  s  +  t]|Es} 


=  exp  [  -  (c^g  +  c2Rg)t]  . 

This  proves  that  the  process  ((Bt,  Rt))  has  jump  function  X  defined 
in  (  7)  above . 


We  next  show  that  for  (i,  j)  e  Eq, 


P  ’  l<Bv+’  ^  "  1»  ^)}  ic.  +  jc, 


and 


P(i,j)KBv+>  <V  =<i>  3  -  !>l 


where 


(BV+’  V>  (BV+t’  RV+t) 


is  the  state  entered  by  the  process  at  the  time  V  of  the  first 
fatality.  It  seems  plausible  that  that  state  must  be  either 
(i  -  1,  j),  if  the  first  fatality  is  suffered  by  Blue,  or  (i,  j  -  1), 
if  the  first  fatality  is  suffered  by  Red. 
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Let 


sb  min ^bsi  ^  *  •••>  BS1  ^ 

and 

SR  =  min|RS^  \  RS^^|  . 

Then 

PU’j)lfflv+,  R^.)  =  (1  -  1,  j)(  =  P(i>j)lv  =  SR| 

=  P<i,j)lsB  >  sRl 

and 

P(1,3)1<V>  V5  ■  (i>  j  -  x)l  =  P<i,;i)tv  =  SB| 

■  <  SR|  • 

Since  SR  is  exponentially  distributed  with  parameter  ic1  and  SR  is 
exponentially  distributed  with  parameter  jc2  and  SR,  SR  are  independent 
(all  these  assertions  hold  with  respect  to  the  probability  measure 
of  course),  it  follows  that 

P<i’j)|SB  =  SRi  =  0 


and  therefore 

P(i,j){(Ev+,  R v+)  ef(i  -  1,  j),  (i,  j  -  1)}|  =  1  . 


We  next  compute  t;hat 


E(i>  j  )[-p 

jc2  / 

0 

2.C1  +  JC 


(i*j)|SR<SB|SR|] 

exp  (  -  jc2v)  exp  (  -  ic  v)  dv 


2 
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where  denotes  expectation  with  respect  to  the  probability  P^’^. 

Similarly, 


P(i’J)  |sB  <  S„|  = 


1C 


1 


R 1 


ic1  +  jc2 


For  each  t  let  denote  the  time  of  the  first  fatal  shot  fired 
after  time  t  .  Then  again  using  the  MmemorylessM  property  of  the 
exponential  distribution,  one  obtains 

P  ,j>Kt+  ’  ^4^  =  ^k’  ^lltl 

=  p  K+>  V)  =  (k>  ^ 


=  P((Bt,  Rt),  (k,  £)), 

where  P  is  defined  by  (8),  for  each  t  and  (k,  i) . 

Hence  ((Bt,  Rt))t>g  has  the  Markov  property,  the  given  jump 
function,  and  the  given  transition  kernel.  That  this  process  then 
has  the  infinitesimal  generator  Q  given  in  (9)  follows  from  the 
theory  of  regular  step  processes,  cf  p.  57. 

Here  are  some  additional  properties  of  this  process  which  are  of 
physical  interest. 

COROLLARY.  The  following  properties  hold. 

1)  With  probability  one,  the  functions  t  -•  B.  and  t  -*  R.  are 
nondecreasing ; 

2)  If  i  /  0,  j^O,  then 

p( 1 » j  )  |  (b  Rt)  =  (0,  0)  for  some  t[  =  0  . 

Property  1)  certainly  makes  sense:  if  there  are  no  reinforce¬ 
ments,  surviving  force  strengths  can  only  decrease.  The  second 
assertion,  which  is  true  even  if  c^  =  c 2  and  initial  force  strengths 
are  numerically  equal,  states  the  nmutual  annihilation  at  time  t  =  <=' 
which  occurs  in  the  deterministic  Lanchester  square  model  is  almost 
surely  impossible  in  this  model. 
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We  wish  to  emphasize  that  although  the  vector  process  ((Bt, 
is  a  Markov  process,  neither  of  the  individual  survivor  processes 

and  is  Markovian,  at  least  with  respect  to  the  family  of 

histories  (Et)t>Q  If  t  is  fixed,  then  the  history  E  contains 
information  other  than  the  current  state  which  is  helpful  in 
predicting  the  future  course  of  the  process  (Bu);  that  is,  in 
predicting  (Bg)  That  information  is,  of  course,  the  current 

number  R  of  Red- survivors .  (It  couldn’t  be  anything  else,  for 
otherwise  ((B  ,  R^))  wouldn’t  have  the  Markov  property.) 

This  model  can  be  extended  to  allow  for  reinforcements  which  arrive 
in  a  particular  manner,  namely,  according  to  compound  Poisson  processes, 
as  we  now  indicate. 

Feller  (1966)  contains  sufficient  background.  For  our  purposes 
a  compound  Poisson  process  is  defined  by  a  positive  number  a  and 
a  probability  distribution  cp  on  N*  =  {  1,  2,  ...}.  Let  Tq  =  0, 

Tp  Tgj  •••  he  the  times  of  arrivals  of  a  simple  Poisson  process 
with  rate  a  and  let  X^,  X2,  ...  be  independent  and  identically 
distributed  as  cp  .  Assume  that  (T  )  and  (X^)  are  independent  and 
define 

M  =  E  X  ,  t  >  0  . 

n:T  <t 

n— 

One  interprets  (Nt)  is  an  arrival  counting  process;  that  is, 
the  number  of  objects  which  have  arrived  up  until  time  t  . 
of  the  times  T^,  T2,  ...  a  random  number  of  objects  arrives, 
distributed  as  cp  . 

Consider  now  the  following  additional  assumptions. 

6.  Reinforcements  to  the  Blue  side  arrive  according  to  a  compound 
Poisson  process  with  parameters  (a-^j  cp^)*  Reinforcements  to  the  Red 
side  arrive  according  to  a  compound  Poisson  process  with  parameters 
(a CP2 ) -  All  reinforcements  immediately  enter  the  combat. 


Nt  is 
At  each 
which  is 
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7.  The  two  reinforcement  processes  are  independent  of  each  other 
and  of  all  firing  processes. 

For  (i,j)  e  E,  let  p^’^)  denote  the  probability  law  of  the  homo¬ 
geneous  square  law  process  with  reinforcements,  governed  by  assumptions 
1-7  above  and  conditioned  on  the  event  {  (B^,  Rq)  =  (i,  j)|. 

(10)  THEOREM.  The  process 

(0,  E,  Et,  (Bt,  Rt),  P(l>;,)) 
is  a  regular  step  process  with 

a)  state  space  (E,  E)  as  defined  for  Theorem  (6); 

b)  jump  function  X  given 

|ic1  +  jc2  +  cl1  +  a2  if  i  >  0  and  j  >  0 

. 

(a-L  +  a2  if  i  =  0  or  j  =  0  ; 


c)  transition  kernel  P  given  for  states  (i,  j)  with  i  >  0  and 
j  >  0  by 


1C- 


P((i,  j);  (i,  j  -  D)  =  - 

P((i,  j);  (i  -  1,  j))  = 


1C1  +  JC2  +  ai  +  a2 

7co 


(12a) 


1C1  +  2  c2  +  ai  +  a2 


a- 


P((i,  j);  (i  +  k,  j))  =  cp1(k)  iCl  ;  jc2  TJ-TV, 
P((i,  j);  (1,  j  +  i)>  ^,(»)  iv  .c*2  a.+a, 

and  for  states  (i,  j)  with  i  =  0  or  j  =  0  by 


a. 


(12b) 


P((i,  j);  (i  +  R,  j))  =  cp-|(k)  - — —  , 

1  +  0-2 


a. 


P((i,  j);  (i,  j  +  JO)  =  cp «(£)  - — - 

A  +  cl. 


k  >  1 


i  >  1 


k  >  1 


l  >  1  ; 


74 


d)  infinitesimal  generator  Q  given  for  states  (i,  j)  with 
i  >  0,  j  >  0  by 

Q((i,  j);  (i,  j  -  1))  =  ic± 

Q((i,  j);  (i  -  1,  j))  =  jc2 

(13a)  Q((i,  j);  (i  +  k,  j))  =  a^Ck) 

Q((i,  j);  (i,  j  +  £))  =  a^>2U) 

Q((i,  j);  (i,  j))  =  -  (ic1  +  jc2  +  a  +  a  ) 

and  for  states  (i,  j)  with  i  =  0  or  j  =  0  by 

Q((i,  j);  (i  +  k,  j))  =  a1cp1(k) 

(  13b)  Q((i,  j);  (i,  j  +  £))  =  a^U) 

Q((i,  j);  (i5  j))  =  -  (ax  +  a2) 

PROOF.  We  employ  the  same  notation  as  in  the  proof  of  the  preceding 
Theorem,  with  the  following  additions.  Let  denote  the  time  of  the 
arrival  of  the  first  group  of  Blue  reinforcements,  let  A_  denote  the 
time  at  which  the  first  contingent  of  Red  reinforcements  arrives .  Let 
Xg,  denote  the  respective  sizes  of  these  two  groups  of  reinforcements. 

The  time 

V  =  inf  jt:  (Bt,  Rt)  ?  (BQ,  RQ)} 

at  which  the  process  first  changes  state  is  now  given  by 

V  =  min|Sg,  S^,  Ag,  A^J 

where  Sg,  are  as  defined  in  the  proof  of  Theorem  (6). 

That  is,  the  first  change  of  state  occurs  in  exactly  one  of  the 
following  ways: 

i)  A  fatal  shot  is  fired  by  Blue  before  Red  fires  a  fatal  shot  and 
before  any  reinforcements  arrive;  that  is 

SB  <  SR5  SB  <  AB’  SB  <  AR 
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ii)  Red  fires  a  fatal  shot  before  Blue  does  and  before  any 
reinforcements  arrive: 

SR  "=  SB>  SR  <  AB’  SR  ■=  SR 

iii)  Some  Blue  reinforcements  arrive  before  any  Red  reinforcements 
do  and  before  any  fatal  shots  are  fired: 

AB  <  SB’  AB  <  SR’  AB  <  AR 

.v)  Some  Red  reinforcements  arrive  before  the  arrival  of  any  Blue 
reinforcements  and  before  any  fatalities: 

AR  <  SB>  AR  <  SR>  *R  <  RB  ' 

By  assumptions  2,  3,  4,  5,  6,  and  7,  the  random  variables  Sg,  Sg, 

A_  ,  are  independent  and  exponentially  distributed  with  means  1/ic  , 

B  R  -L 

l/jc0,  l/a1  and  l/a^  3  respectively.  Therefore  V  is  exponentially 
distributed  with  mean  l/(ic  +  jc 2  +  +  a2)  =  ± , j ) ;  in  other  words 

(  14)  P^1*^ ^  ](Bu,  Ru)  =  (i,  j )  for  0  <  u  <  t}  =  exp[  -  X(i,  j )t] 

for  all  (i,  j)  and  t  ,  where  \  is  given  by  (11). 

Moreover,  since 

P<1,;i)!SB  =  SR  °r  SB  =  AB  °r  SB  =  AR  °r  SR  =  AB  °r  SR  =  AR  0r  AB  =  ARI  =  0 
it  follows  that 

p(1’j)kBv+’  V’  e  = 1 

where 

Eij  =  i(i  -  1,  j ),  (i,  j  -  1)  \  U  |  (k,  j):  k  >  i} 

U  'f(i,  A):  l  >  jj. 

That  is,  two  or  more  simultaneous  fatalities  (on  either  the  same  or 
opposite  sides)  or  a  simultaneous  fatality  and  arrival  of  reinforcements, 
fail  to  occur  almost  surely. 
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Proceeding  as  in  the  proof  of  Theorem  ( 6)  we  see  that 

p(1J)l<Bv+’  v*  =  (i>  i  -  1)l 

'j)|sB  <  SR,  SB  <  Ag,  SB  <  Br) 

_  _ _ 

iGl  +  ~*C2  +  ai  +  a2  ’ 

that 

■  p(i,:,)  l<V»  RV+)  =  (i  “  1> 

p^  ’^|sR  <  SB,  SR  <  Afi,  SR  <  Ar| 
jC2 

ic  ■t  jCj  ^  ct  ^  01 2 

that 

l<Bv+»  V5  et(k’  j):  k  >  iil 

=  p(iJ)K  <  SB>  AB  <  SR-  AB  <  ArI 

=  al 

ic1  +  j^2  +  oc 2 

and,  finally,  that 

P<1,;|)I(BV+,  Rv+)  e|(i,  iy.  je  >  j n 

=  p(iJ,tAR  <  SB’  ar<  sr>  ar  <  ab( 
a  2 

"  Ic7  +  jc2  +  a±  +  a2  ’ 

Let  us  consider  the  latter  two  expressions  in  more  detail.  Since 
(B^+,  Ry  )  =  (k  +  i,  j)  with  k  >  0  if  and  only  if  the  first  event 
causing  a  change  of  state,  is  the  arrival  of  a  group  of  Blue  reinforcements 
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which  number  exactly  k,  we  have 


p<i,jVv+’  v5  =  <i  +  k> 

_P'  ’^®B  <  SB’  RB  <  SR’  <  RR’  *B  =  * 

By  properties  of  the  compound  Poisson  process  of  Blue  reinforcements 
and  assumption  7,  is,  with  respect  to  the  probability  P^’-^, 
independent  of  A  ,  S_,  SD,  and  A_.,and  therefore 

D  D  K  K 

p(1’j)|AB<SB,  Rb<Sr>  Ab<Sr>  XB=k| 

Cl5a)  =  pCijl  )|^b  <  sbj  ^  <  Sr,  Ab  <  SR}  P(l,j){xB  =  k( 

al 

_  icl  +  Jc2  +  oo-l  +  a2  • 

Through  entirely  analogous  reasoning  we  are  able  to  obtain 

p(i,;i)l<v>  V5  =  (i>  j  + 

(15b)  a9 

=  , - A -  cp9(k) 

1C1  +  3°2  +  ai  +  a2  1 


for  each  z  >  0. 

We  may,  on  the  basis  of  (12),  summarize  (15)  in  the  form 
(16)  P(i’j)l(Bv+,  Ry+)  =  (k,  A) |  =  P((i,j),  (k,  A)) 

for  all  (i,  j),  (k,  Z)  e  E. 
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Using  standard  arguments  concerning  the  exponential  distribution, 
one  extends  (14)  to 

P(1,;|)|(BU,  Ru)  =  (Bs,  Rs)  for  s  <  u  <  s  +  t||s| 

=  exp  r  -  \(BS,  Rs)t] 

and  (16)  to 

pC i > j  ) j  Ry  )  =  (k, 

s  s 

=  P((BS,  Rg);  (k, 

where 

\?s  =  infju  >  sr  (Bu,  Ru)  1  (Bg,  Rg)} 

is  the  time  of  the  first  change  of  state  of  the  process  after  time  t. 

These  last  two  statements  combine  to  complete  the  proof  of  the 
Theorem. 

Note  that  the  process  with  reinforcements  has  no  absorbing 
states . 

Other  types  of  assumptions  concerning  reinforcements  are  possible, 
but  would  be  more  difficult  to  handle  and  are  almost  certain  to  destroy 
the  Markovian  property  of  the  process,  although  if  the  resultant  model 
were  a  semi-Markov  process  it  might  still  be  tractable. 


*>lEsl 

*>>. 
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S2  Heterogeneous  Square  Law  Process 


For  heterogeneous  processes  we  need  to  construct  a  new  sample 
space.  Throughout  we  assume  that  there  are  M  Blue  weapon  types 
(the  well-established  terminology  "weapon  types",  as  opposed  to  "types 
of  combatants",  will  be  used  in  all  our  discussions  of  heterogeneous 
processes)  and  N  Red  weapon  types.  There  is  no  assumption  that  at 
least  one  of  each  type  of  weapon  be  present. 

We  take  as  our  state  space,  therefore,  the  set 

E=NxNx...xN  (M+N  times) 

together  with  the  discrete  a-algebra  E. 

As  sample  space  fi  we  again  have  the  family  of  all  mappings 
oj  =  (cup  ufy,  u/p  iDjJ)  of  mappings  of  [0,  ®)  into  0  which 

are  right  continuous  and  have  left-hand  limits  everywhere.  We  define 
coordinate  vector-valued  random  processes  (gt)t>Q  and  (E1_)t>  q  by 

B^_(uj)  —  (uj-^(t),  ...,  U^j(t)) 

Rt(w)  =  (U){(t),  ...,  (JUjJ(t))  , 

and  coordinate  processes  (B^)^,  .  . . ,  (B™)t>Q,  ^Rt^t>0’ ’  ‘  ‘ ’^Rt^t>0  by 

BR(uj)  =  u^Ct)  ,  1  <  k  <  M 

R*(w)  =  uu'(t)  ,  1  <  l  <  N 

so  that,  in  particular, 

B.J-  =  (B^_,  •  •  •  j  B^ ) 

and 

=  •  •  • »  ) 
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for  each  t  .  B.  is  the  number  of  Blue  type  k  weapons  surviving  at 

^  a 

time  t,  Bt  the  Blue  force  surviving  at  time  t,  the  number  of  type  l 
Red  weapons  surviving  at  time  t,  and  Rt  the  surviving  Red  force  at 
time  t. 

As  before,  we  put 

it  =  ctK§s>  -s):  0  <  s  <  t\ 

and 

E  =  aj(Bs,  Rs):  s  >  0|  . 

£  is  the  history  of  the  attrition  process  up  until  time  t  and  F  its 
entire  history 

M  N 

We  denote  points  in  E  in  the  form  (x,y)  where  x  e  N  and  y  e  N  ; 
x  thus  is  a  state  which  the  vector  of  Blue  surviving  weapons  may 
assume  and  similarly  for  y  and  Red.  x^  corresponds  to  surviving 
Blue  type  i  weapons  and  y^  to  surviving  Red  type  j  weapons . 

For  each  (x,y)  e  E,  let  P^x,y^  denote  the  probability  law  on 
(03  E)  of  the  attrition  process  governed  by  assumptions  1-6  of  the 
Process  S2,  when  (Bq,  Rq)  =  (x,  y). 

(17)  THEOREM.  Under  the  S2  assumptions, 

(0,  E,  Et,  (Bt,  Rt),  P(x’y)) 


is  a  regular  step  process  with 

a)  state  space  (E,  §)  as  defined  above; 

b)  jump  function  \  given  by 


Ux,y) 


.  E  Z  [c  (i, j )y  +  c  ( j , i )x.]  if  x  ^  0,  y  ^  0 

i:x.>0  j:y->0  J 

i 


0 


if  x  =  0  or  y  =  0 
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c)  transition  kernel  P  defined  for  states  (x,y)  with  x  ^  0  and 
y  /  0  by 

M 

E  c1(j,i)xi 

P((x,y);  (x;yu...,y^l . y„»  =  —\<x,y) - 

for  j  =  1,  N. such  that  y_.  >  0  and  by 

N 

E  c2(i,j)yj 

PCC^>y)j  C  ...jX^-ij...j  x^jj  y ) )  _  ~  \^x  y) 

for  all  i  with  1  <  i  <  M  and  x^  >  0; 

d)  infinitesimal  generator  Q  given  for  states  (x,y)  with  x  ^  0 
and  y  /  0  by 

M 

!  Z  c  (j,i)x.  if  y  >  0 
]i=l  3 

QC ( x,  y )  j  —  j 

I  0  if  y.  =  0 

l  3 

Q((x,y);(x,y))  =  -  Z  Z  [c  (i,j)y.  +  c  (j,i)x.] 

i:x^>0  j:yj>0  3 

N 

I  E  c?(i,j)y.  if  x.  >  0 
li=l  3 

Q((x,y);.(x1,...,xi-l,...,xM;y))  =  J 

10  if  x.  =  0 

I  l 

and  for  states  (x,y)  with  x  =  0  or  y  =  0  by  Q((x,y);  £)  =  0  for  all 

C  e  E* 
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PROOF.  We  recall  that  for  q  =  1,2  and  appropriate  k,£  we  have  defined 


cq(k,£)  =  rq(k,l)pq(k,jD  . 

For  i  =  1,  _ ,  M  let 

SR(i)  =  infjt:  B^  /  Bg  f 

be  the  time  at  which  the  first  fatality  to  a  Blue  type  i  weapon  occurs. 

If  there  is  no  such  time,  and  in  particular  if  Bg  =  0,  we  put  SR(i)  =  +  ®. 
Similarly,  let 

SB(j)  =  infjt:  Rj  ±  R^}  ,  j  =  1,  ...,  N 

be  the  time  of  the  first  fatality,  if  there  is  one,  to  a  Red  type  j 
weapon. 

We  claim  that 

P(x’y)jsR(i)  >  t,  i  =  1,  ...,  M,  SB(j)  >  t,  j  =  1,  ...,  N( 

(18) 

=  exp  (  -  t  E  E  [c«(i, j )y .  +  c  (j,i)x,]) 

i:x^>0  j:yj>0  ^ 

from  which  it  follows  that,  if  we  define 

V  =  infjt:  (Bt,  Rt)  *  (Bg,  Rg)| 

then  since 

V  =  minjSg(l),  •••}  Sg(N),  SR(1),  •••>  SR(M)j  > 

one  has 

p(x,y)j V  >  t}  =  exp  [  -  X(x,y)t] 
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Now,  how  are  we  to  prove  (18)?  If  all  kills  were  hypothetical 
only,  then  for  a  given  j  and  i  ,  the  times  of  hypothetical”  kills 
of  type  j  Red  weapons  by  type  i  Blue  weapons  would,  by  assumptions  2 
and  4  form  a  Poisson  process  with  rate  c^(j,i)x^,  so  long  as  y.  >0. 

Hence  by  hypothesis  6,  "kills”  of  type  j  Red  weapons  occur  hypothetically 
at  the  arrival  times  of  a  Poisson  process  with  rate 

Z  c1(j,i)xi  ; 
i 

in  particular  the  time  s£(j)  of  the  first  such  hypothetical  kill  is 

a  -1 

exponentially  distributed  with  mean  (  E  c^(j,i)x. )  if  y.  >0  and 

i  1  J 

mean  +  ®  if  y.  =0.  Moreover,  by  hypothesis  3,  S?(l),  ...,  S*(N)  are 

3  (x  v)  a  B 

independent  with  respect  to  P  ’  . 

If  S£(i)  similarly  stands  for  the  time  of  the  first  hypothetical 
kill  of  a  Blue  type  i  weapon,  then  S^(i)  is  exponentially  distributed, 

with  respect  to  the  probability  P^X’^,  with  mean 

(  (E  c  (i,j )y. )-1  if  x.  >  0 

j  3  *  J  1 

J  +  oo  if  x^  =  0 

Assumption  3  assures  us  that  S|j(l),  ...,  S||(M)  are  independent  and 
assumption  6  implies  that  |s^(l),  ...,  S*(N)j  and  |S|f(l),  ...,  S|;(M)} 
are  independent. 

Furthermore,  even  though  by  definition 

V  =  min|SB(l),  ...,  Sb(N),  Sr(1),  SR(M)f 

it  is  also  true  that 

V  =  minjs*(l),  S*(N),  S*(l),  S*(M)} 

since  the  times  of  the  first  actual  and  hypothetical  kills  necessarily 
coincide. 
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Thus 


p(x,y)|v  >  t\  =  P(x’y){s*(l)  >  t,  ...,  S*(N)  >  t  , 

S*(l)  >  t,  S*(M)  >  t} 

=  P(X,y)  {S*(l)  >  t}  ...  P(X,y){S*(N)  >  t} 

x  P(X’y)js*(l)  >  t|  ...  P(X,y)|s*(M)  >  tf 

=  exp  [  -  t  E  E  c.(j,i)x.] 
j:yj>°  1 

x  exp  [  -  t  EE  c  (i,j)y.] 
i : x^>0  j  J 

which  completes  the  proof  of  (18). 

We  note  that  if  j  ^  l,  i  ^  k 

P(X,y)  (S*(j)  =  S*(£),  S*(j)  <  S*(£)  <  <0}  =  0 

and 

P(X,V)  |S*(i)  =  S*(k) ,  S*(i)  <  co,  s*(£)  <  »}  =  0 

and  that 

P(X,y)  |s*(j  )  =  S*(i),  S*(j  )  <  s*(i)  <  “I  =  0 

for  all  i,j,  by  independence  and  exponential  distributions  of  the 
random  variables  involved. 

Since 

p(  ,y)|(By+,  Ry^_)  =  (x;y1}  •••»  yj-^’  ••*’  yN^ 

=  P(x’y)js*(j)  <  S*(£),  i  t  j;  S*(j)  <  S*(i)  for  all  i} 
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it  follows  that 


p(x>v)  I (By+,  Ry+)  e|(w,z):  wi  =  x^-1  for  exactly  one  i  and 

otherwise  w  =  x  and  z  =  y  or  z.  =  y.-l  for 
exactly  one  j  and  otherwise,  z  =  y  and  w  =  x} 

Again  by  independence  and  exponential  distributions 

P(X’y)  !s*(j)  <  S*u),  &  *  j;  S*(j)  <  S*(i)  for  all  i} 

0  if  y.  =  0 

3 

=  Z  c1(j,i)xi 

- - -  if  y.  >0 

Z  E  c  (j,i)x.  +  E  E  c  (i,j)x  J 

j:y.>0  i  i:x.>0  j  J 

■  j  i 

and  similarly, 

p(x,y){(§v+’  ^  =  (xi’  •••>  xi-1»  •••>  v  y)l 

=  P(x’y)  |s*(i)  <  S*(k),  k  /i;  S|(i)  <  S*(j)  for  all  j[ 

0  if  x.  =  0 

i 

=  z  c2(i, j )yj 

— - - ^ -  if  x.  >  0 

Z  E  c  (j,i)x  +  Z  Z  c  (i,j)y. 

j:y.>0  i  i:x.>0  j  J 

3  i 

Hence  we  have  shown  that 

p(X’y)  |CBV+,  R v+)  =  ?|  =  PC(x,y);  ?) 

for  all  ^  e  E,  where  P  is  the  transition  kernel  defined  in  the 
statement  of  the  Theorem.  [1 
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S3a  Heterogeneous  Square  Law  Process 

We  continue  to  assume  that  there  are  M  Blue  weapon  types  and 
N  Red  weapon  types,  and  that  the  canonical  stochastic  process 

(0,  E,  Et>  (Bt,  Rt)) 

f  x  y  ) 

is  as  described  in  the  preceding  section.  PK  will  now  denote  the 

probability  law  of  the  attrition  process  governed  by  the  family  of 
five  assumptions  for  the  Process  S3a,  subject  to  the  initial  conditions 


.M 


B.q  =  x  e  N 

N 

R0  =  y  e  If  . 


(19)  THEOREM.  Under  the  assumptions  of  Process  S3a,  the  collection 

(n,  E>  Et,  (Bt,  Rt),  p(x>y)) 


is  a  regular  step  process  with 

a)  state  space  (E,  E)  as  defined  in  the  description  of  Process  S2; 

b)  jump  function  X  given  by 


M  N 

I  Z  .Z  [x.p,(i,j)r  (j)y.  +  y .p, ( j ,i)r1(i)x.]  if  x  #  0,  y  /  0 
i=l  j=l  1  3  3 

\C,x,y)  =( 


0  if  x  =  0  or  y  =  0 


.  where 


and 


M 

x  •  =  x .  /  E  x.  ,  i  =  l,  . . . ,  M 
1  k=l 

N 

yH  =  yJ.Z  y,  >  j  =  1,  ...»  n  ; 

3  3  t=i  1 
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these  quantities  are  defined  only  if  the  denominator  is  positive; 

c)  transition  kernel  P  given  for  states  (x,y)  e  E  such  that 
x  ^  0  and  y  ^  0  by 


p((*,y);  (x;y15 


,  yj  -  !,  •••>  yN)) 

M 


E  P-,(j,i)r  ,(i)x. 
i=l  _ 

l(x,y) 


j  =  1,  . . . ,  N 


pc(x5y)i  ( x^,  • • • 5  xi_ • • • 5  y)) 


N 


=  x. 


E  P2(i,j)r2(j)y, 
3=1 


i  =  1  M  • 


i  \(x,y) 

d)  infinitesimal  generator  Q  given  for  states  (x,y)  with  x  ^  0 


and  y  /  0  by 


M 


Q((x,y),  (x;y15  ...,  y_.-l,  ...,  yN))  =  y_j  E_^  P-^ j  ,i)r1(i)x± 

M  N 

Q((x,y),  (x,y) )  =  -  E  E  [x.p9(i,j)r  (j)y  +  y.p  (j,i)r ,(i)x.] 

i=l  j=l  J  J 

N 

Q((x,y),(x^,  •••)  x^-lj  •  ••>  xM;y))  —  x^  E  P2C^->3  * 


If  x  =  0  or  y  =  0,  Q((x,y),  a)  =  0  for  all  a  e  E. 

PROOF.  As  usual  consider  the  hypothetical  situation  in  which  no 
weapons  are  actually  killed,  but  only  "potential"  kills  are  recorded, 
in  terms  of  the  times  at  which  they  occur.  Following  the  notation  of 
the  preceding  proof  we  let  S*(j)  denote  the  time  of  the  first  hypo- 
thetical  kill  of  a  Red  type  j  weapon,  j  =  1,  . . . ,  N  and  for  i  =  1, 
...,  M  we  denote  by  S|;(i)  the  time  at  which  the  first  hypothetical 
kill  occurs  of  a  Blue  type  i  weapon.  Then  the  time 
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V*  =  min|Sg( 1),  ...,  S*(N),  S*(l),  S*(M)f 

of  the  first  hypothetical  kill  coincides  with  the  time 

V  =  inf  jt:  (Bt,  Rfc)  t  (BQ,  RQ)f 

of  the  first  actual  kill,  which  is  also  the  time  at  which  the  attrition 
process  first  changes  state. 

For  i  =  1,  M  and  j  =  1,  N  let  K|f(i,j)  be  the  time  of  the 

first  hypothetical  shot  which  is  fatal  to  a  Blue  type  i  weapon  and  was 
fired  by  a  Red  type  j  weapon.  Similarly,  let  K*(j,i)  be  the  first 
time  at  which  a  Red  type  j  weapon  is  hypothetically  killed  by  a  Blue 
type  i  weapon.  Then 

S*(i)  =  min{K*(i,l),  ...,  K*(i,N)f 

and 

S*(j)  =  min{K*(j,l),  ...,  Kj(j,M)( 

for  all  appropriate  i  and  j  . 

Consider  now  a  fixed  Blue  weapon  of  type  i.  The  times  of  hypo¬ 
thetical  shots  which  it  fires  form  a  Poisson  process  with  rate  r1(i). 
Further,  by  assumption  3 

p(x>y)|a  hypothetical  shot  fired  by  a  Blue  type  i  weapon 
is  directed  at  some  Red  weapon  of  type  j } 

(2°)  N  • 

=  h  =  V  £  y.  > 

this  is  also  the  probability  that  an  actual  shot  fired  at  a  Red  force 
of  composition  y  is  directed  at  a  target  of  type  j .  A  similar  ex¬ 
pression  holds  for  shots  and  hypothetical  shots  fired  by  Red  weapons. 

Hence  hypothetical  shots  fired  by  a  particular  Blue  type  i  which  are 

directed  at  Red  weapons  of  type  j  form  a  Poisson  process  with  rate  r, (i)y. 

■*-  .  D 
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with  respect,  of  course,  to  the  probability  P^x,y\  which  is  assumed 
fixed  for  the  purposes  of  this  discussion.  The  fatal  shots  fired  by 
this  weapon  at  Red  type  j  weapons  thus  constitute  a  Poisson  process 
with  rate  r^(i )y^p^( j , i ) .  By  assumption  6,  therefore,  hypothetically 
fatal  shots  fired  at  Red  type  j  weapons  by  all  type  i  Blue  weapons 
form  a  Poisson  process  with  rate  x^r^(i)y_.p^( j ,i). 

Another  application  of  assumption  6  shows  that  the  times  of  all 
hypothetically  fatal  shots  fired  at  Red  type  j  weapons  form  a  Poisson 
process  with  rate 


E  ><iriCi)yiP1C3 ,i) 
i  J 


and  therefore 

(21)  P(x,y)is*(j)  >  t(  =  exp[  -  ty^  S  P1(j,i)r1(i)xi]  . 
Similarly,  since 

p(X,y)|K*(i,  j  )  >  t}  =  exp[  -  t^PjCi,  j  )r2(  j  )y_.] 

for  each  i  and  j  yet  another  application  of  the  independence 
assumption  6  yields 

(22)  p(x>y)js*(i)  >  t}  =  exp  [  -  tx.  E  p2(i, j )r2( j )y . ]  . 

j 

From  these  two  expressions 

p(*>y)|v*  >  t|  =  exp[  -  X(x,y)t] 
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follows  and  hence 

pC x? .V ) j V  >  t}  =  exp[  -  X(x,y)t]  . 

Proceeding  exactly  as  in  the  proof  of  Theorem  (17)  we  see  that 
^  5  I ( Sy+  5  By+)  -  (x;y^,  •••,  yj~l>  •••5 

=  p(x’y)|v*  =  S*(j)| 

(23) 

=  P(X’y)|s*(j)  <  S*U)  for  i  I  j,  S*(j )  <  S*(i)  for  all  ij 

yj  E  P1(j,i)r1(i)xi 
X(x,y) 

and,  analogously,  that 

p(X’y)|(Bv+,  Ry+)  =  (x-L,  x.-l,  ...,  xM;y)f 

=  P(x>y)j v*  =  S*( i ) \ 

(24) 

=  P(x’y)|s*(i)  <  S*(k)  for  k  +  i,  S*(i)  <  S*(  j  )  for  all  j} 
xi  ? 

=  x^Ty)  '  LI 


93 


S3  Heterogeneous  Square  Law  Process 

The  notation  is  exactly  that  of  the  previous  two  processes. 
pO,y)  now  denotes  the  probability  law  of  the  attrition  process 
governed  by  the  family  of  assumptions  of  Process  S3,  subject  to  the 
initial  conditions  (BQ,  RQ)  =  (x,y). 

(25)  THEOREM.  Under  the  assumptions  for  Process  S3,  the  collection 

(0,  E,  Et,  (Bt,  Rt),  P(X’y)) 


is  a  regular  step  process  with 

a)  state  space  (E,  |)  as  defined  above; 

b)  jump  function  \  given  by 

M  N 

X(x,y)  =  E  E  [Yi(y,j)p-,(j,i)r1(i)x.  +  q-(x,i)p?(i,j)r5(j)y  ] 
i=l  j=l  x  ±13  ^  z  3 

(we  assume  that  y^  =  0  implies  Y^(y,j)  =  0  for  all  i); 

c)  transition  kernel  P  given  for  states  (x,y)  with  x  ^  0  and 
y  /  0  by 


P((x,y);(x;y1, 


V1' 


M 

S  Y±(y, j)p1(j,i)r1(i)xi 
•••’  yN)}  =  X(x,y) 


P((x,y);(x1} 


. ..,  x±-l, 


•  •  •  9 


vy)) 


E  r|.(x,i)p2(i,  j  )r2(j  )y  . 
3  J  J 

- TU^y) - 


d)  infinitesimal  generator  Q  given  for  states  (x,y)  with  x  /  0 
and  y  ^  0  by 


Q((x,y);(x;y1} 


yj-i,  •••>  yK» 


M 

S  Y,(y,j)p1(j,i)r,(i)x. 
i=l  1  x 
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M  N 

Q((x,y);(x,y))  =  -  E  E  [Y .(y, j )p,( j,i)r  (i)x.  + 

i=l  j=l  1  1  11 

+  rij(x,i)p2(i, j  )r2(j)y_.] 


Q((x,y);(x1, 


N 


x .  “1, 


^jy))  =  E  Ti.(x,i)p2(i,  j)r  2(j)y. 

j=l  J  J 


If  x  =  0  or  y  =  0,  Q((x,y),a)  =  0  for  all  a  e  E. 

PROOF.  The  proof  of  Theorem  (19)  requires  only  notational  changes. 
By  the  new  assumption  3,  (20)  becomes 


(20*)  p(x>y)  I  a  hypothetical  shot  fired  by  a  Blue  type  i  weapon 
is  directed  at  some  Red  weapon  of  type  j  } 

=  ^(yj) 

so  that  hypothetical  shots  fired  by  a  given  Blue  type  i  weapon  at 
Red  type  j  weapons  form  a  Poisson  process  with  rate  r1(i)Yi(y, j ) 
and  hypothetical  fatal  shots  fired  by  this  type  i  weapon  at 
Red  type  j  weapons  form  a  Poisson  process  with  rate 
ri(i)Yi(y, j )pp( j ,i),  in  which  case  for  (21)  one  obtains 


(21*)  P(X,y)  |  S*(j)  >  t } 

=  exp[  -  t  E  Yi(y,j)p1(j,i)r1(i)xi]  ; 
i 

analogously,  (22)  becomes 


(22*)  P(X’y)  |  S*(i)  >  t| 

=  exp  [-  t  E  n.(x,i)p2(i,j)r2(j)y.]  . 

j 
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In  the  same  manner  one  obtains  as  analogs  of  (23)  and  (24)  the 
following: 


(23*) 

p(x,y){s*(j)  <  s*(£)  for  l  +  j,  S*(j)  <  S*  (i)  for  all  if 

E  7i(y, j)p1(j,i)r1(i)xi 

- Wtt - 

for  j 

=  1,  . . . ,  N  and 

(24*) 

P(X,y)|s*(i)  <  S*(k)  for  k  *  i,  S*(i)  <  S*(j)  for  all  jf 

E  ri.(x,i)p2(i,j)r  (j)y 

3 

X(x,y) 

respectively. 

The  Theorem  now  follows.  If 

EXAMPLES.  The  fire  allocation  scheme  of  this  process  is  quite 
general,  but  also  entirely  abstract.  Here  are  three  specific  cases. 


1. 

Uniform  Fire  Allocation  occurs  when 

and 

N 

Y.(y, j)  =  y-  =  y •/  s  yt 

1  3  3  z=i  1 

M 

cp.(x,i)  =  x.  =  x  /  E  x, 

3  k=l 

for  all  i  and  j  and  in  this  case  the  Process  S3  reduces  to  the 
process  S3a.  It  becomes  clear,  therefore,  that  Process  S3a  is  a 
Lanchester  square  process. 

2.  Priority  Fire  Allocation  can  occur  in  the  following  manner. 
Consider  a  fixed  type  of  Blue  weapon,  say  type  i,  and  suppose  that 


97 


the  Red  weapon  type  1,  . ..,  N  are  ranked  in  some  order  j.  ,  . ..,  j 
(here  { .  ..,  j  }  is  simply  a  permutation  of  the  set  |l,  . ..,  N|) 
in  such  a  manner  that  a  Red  weapon  of  type  j  will  never  be  fired 

Xj 

upon  when  a  weapon  of  type  j  ^  is  present.  One  models  this  by 
putting 


^  y  >  j  £*)  ^ 


where 


l*  =  infjjfc:  y.  >  0|  . 

3  i 

Other  more  flexible  and  possibly  more  realistic  priority  schemes  can 
certainly  be  developed. 

3.  One  may  also  derive  fire  allocations  from  physically  meaningful 
axioms.  For  example,  consider  the  following  hypotheses  on  the  family 
|y^|of  fire  allocation  distributions 

i)  Yi(y,j)  =  0  if  and  only  if  y^  =0; 
ii)  For  all  y  and  z 

N  y 

Y,(y  +  z,j)  =  Y.(y,j)  [  E  ■  j  -  VAy  +  z)U)l 

1  1  t=i  zi  1 

N  z 

+  VA  Z,j)  [  E  j—  Y.Cy  +  Z)U)1 

1  1=1  yi  +  zl  1 

for  all  j  =  1,  . . . ,  N. 

The  interpretations  of  all  these  assumptions  are  physicallly 
meaningful.  The  first  states  that  a  weapon  type  not  present 
receives  no  fire  but  that  every  weapon  type  present  receives 


a  positive  fraction  of  the  fire.  The  second  hypothesis  is  best 
explained  step-by-step.  Consider  the  effect  of  combining  two 
target  forces  y  and  z  into  the  single  force  w  =  y  +  z. 

Y^(y  +  z,  JL)  is  the  proportion  of  fire  directed  at  the  combined 
force  w  which  is  allocated  to  type  i  weapons.  If  this  is 
further  allocated  among  the  type  SL  weapons  from  the  two  com¬ 
ponent  forces  (y  and  z)  in  proportion  to  the  relative  numbers 
of  such  weapons  present  then 


y*  +  Z* 


*  (y  +  z,  i) 


is  the  fraction  of  fire  directed  at  the  combined  force  that  is  allo¬ 
cated  to  type  l  weapons  originally  part  of  the  y- force.  Thus 


N 


SL 


ay  =  E  v  +  z 
y  1=1  y l  i 


T  (y  +  z,  jQ 


is  the  proportion  of  fire  allocated  to  weapons  of  types  originally 
and  the  y-force  and 


N 


SL 


CL-  -  E  - 
z  1=1  yl  +  zi 


Y  (y  +  Z,  SL) 


the  fire  allocated  to  weapons  originally  in  the  z-force.  But  the 

fire  represented  by  should  be  allocated  among  weapon  types 

according  to  the  distribution  Y  (y,  •)  and  similarly  for  a  so  that 

z 

one  should  have 

Y  (y  +  z,  •)  =  ayY  (y,  •)  +  azY  (z,  •) 

which  is  exactly  the  second  assumption.  If  an  assumption  like  this 
did  not  hold,  consistency  problems  would  arise,  with  fire  allocation 
dependent  on  names  given  targets  rather  than  only  numbers  of  targets. 


99 


(For  example,  arbitrarily  splitting  a  class  of  n  indistinguishable 
weapons  into  two  subclasses  of  n1  ’’Type  A”  and  n  -  n^  ’’Type  B” 
weapons  would  change  the  fire  allocation). 

We  have  conjectured,  and  J.  Blankenship  (1973)  has  proved,  using 
the  notion  of  invariant  measures  for  Markov  matrices,  that  any  family 
{y  |  of  allocation  distributions  satisfying  i),  ii),  iii)  above  is 
necessarily  of  the  form 


a  .y . 

Y(y,j)  =  ■  -N-  -  -3 . 

z  aiyi 
1=1  1  L 


where 


=  Y( 1, j )  ,  j  =  1,  ...,  N 

and  1  =  (l,  ...,  l).  It  is  possible,  moreover,  to  express  the  Y 
and  terms  of  Y  (y,  *)  for  any  fixed  y  such  that  y^  >  0  for  all  j. 
Thus  if  one  accepts  the  hypotheses  i),  ii),  iii)  then  the  fire  alloca¬ 
tion  for  every  target  force  can  be  computed  from  that  for  a  ’’base” 
force. 

We  remark  in  conclusion  that  even  if  different  weapons  may  employ 
qualitatively  different  means  of  allocations,  the  model  remains  valid. 
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LI  Homogeneous  Linear  Law  Process 


We  return  to  the  sample  space  and  canonical  stochastic  process 
previously  established  for  the  homogeneous  case.  Thus,  the  state 
space  is 


E  =  N  x  N 

together  with  the  discrete  a-algebra  E,  Q  is  the  set  of  mappings 
uu  =  (uu^>  UU2 )  from  [0,  »)  into  E  which  are  right-continuous  with 
left-hand  limits  everywhere, 

B^Ccu)  =  u)^(t)  e  E 

Rt(u))  =  (u2(t)  e  E 

for  each  t  >  0, 


Et  =  °((BU,  Ru)  •*  0  <  u  <  t) 

is  the  history  of  the  attrition  process  up  until  time  t  and 

F  =  a((Bu,  Ru)  :  u  >  0) 


the  entire  history. 

From  each  (i,j)  e  E  we  denote. by  P^1’3  ,  as  usual,  the  probability  of 
the  attrition  process  conditioned  on  the  event  {(Bg,  Rg)  =  (i,j)} 
and  this  time  governed  by  the  five  assumptions  listed  in  the  main 
text  for  the  Process  LI. . 

(26)  THEOREM.  Under  the  assumptions  of  Process  LI,  the  collection 

(n,  E,  Et>  (Bt,  Rt),  P(i’j)) 
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is  a  regular  step  process  with 

2 

a)  state  space  E  =  N  ; 

b)  jump  function  \  given  by 

=  ij(kx  +  k2); 

c)  transition  kernel  P  given  for  states  (i,j)  with  i  >  0  and 
j  >  0  by 

k. 


P((i,j),  (i,j  -  1))  = 


k1  +  jq 


P((i,j),  (i  -  l,j))  =  r- — — t —  ; 

K1  +  k2 

d)  infinitesimal  generator  Q  given  by 
Q((i,j )>  (i, j  -  1))  =  kjij 
Q((i,j),  (i> j ))  =  -  (kx  +  k2)ij 
Q((i,j),  (i  -  1,  j))  =  k2ij  . 

We  remind  that  for  l  =  1,  2 

k„  =  s„q„ 
j l  SL  SL 

where  1/s^  is  the  mean  time  required  for  a  particular  Blue  weapon  to 
detect  a  particular  Red  weapon  and  q^  the  conditional  probability  that 
a  Blue  weapon  destroys  a  Red  weapon  (in  a  one-on-one  engagement,  the 
only  sort  of  engagement  permitted  in  this  model,  which  is  also  assumed 
to  occur  instantaneously)  given  detection  and  attack.  s2  and  q2  are 
analogous  descriptions  of  Red  weapons. 

PROOF.  We  again  use  the  device  of  hypothetical  kills,  which  is 
useful  because  it  preserves  independence  and  stationarity  properties 
which  are  not  valid  for  the  actual  attrition  process,  yet  allows  us 
to  make  statements  about  the  actual  process  because  the  time  of  the 
first  hypothetical  kill  is  also  the  time  of  the  first  actual  kill 
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(namely,  the  first  time  at  which  the  attrition  process  changes  state) 
and  the  weapon  killed  at  the  first  hypothetical  kill  is  also  the  first 
weapon  actually  killed.  Thereafter  the  two  processes  diverge  from 
one  other  and  the  hypothetical  process  is  no  longer  a  useful  means 
of  dealing  with  the  actual  process. 

We  therefore  denote  by  T*(k,£)  the  time  at  which  the  k  Blue 

weapon  (first)  hypothetically  detects  and  kills  the  l u  Red  weapon 

(not  every  detection  results  in  a  kill  unless  q,  =  1)  and  let  . 

til 

T£(£,k)  be  the  time  at  which  the  £  Red  weapon  first  hypothetically 
K  • 

kills  the  k  Blue  weapon.'  Then  according  to  assumptions  2,  3,  4,  5 
and  the  Random  Sampling  Theorem  for  Poisson  processes,  the  random 
variables  T*(l,l),  T*(l,2),  ...,  T*(l,j),  T*(2,l),  T*(2,2),  ..., 
T*(2,j),  ...,  T*(i,l),  T*(i,2),  ...,  T*(i,j)  are  independent  and 
identically  exponentially  distributed  with  mean  1/k^.  Symmetrically, 
Tr(  1,1),  •••»  Tg(l,i),  •••»  T*(j,l),  ...,  T*(j,i)  are  independent 
and  identically  exponentially  distributed  with  mean  l/kg*  Further¬ 
more,  these  two  families  of  random  variables  are  independent.  All 
of  these  statements  hold  with  respect  to  the  probability 

Also  with  respect  to  P^’^\  the  time 

V*  =  min  min{T*(k,£),  T*(£,k)} 
k<i  B  R 

£<j 

of  the  first  hypothetical  fatality  and  the  time 

V  =  inf |t:  (Bt,  Rt)  ¥■  (BQ,  RQ)} 

of  the  first  actual  fatality,  coincide  almost  surely. 
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By  independence,  therefore, 


p(1>j)jv*  >  t} 


for  1  <  k  <  i,  i  <  £  <  j} 
=  (k,4)  >  t  for  1  <  k  <  i,  1  <  A  <  j} 

x  p( i >  j ) | £ } k )  >  t  for  1  <  k  <  i ,  1  <  £  <  j } 

=  (P(1>J){T*  (1,1)  >  t}  )1J  x  (P(l,j)|T*(l,l)  >  tf)1D 

=  exp  [  -  t(kx  +  k2)ij]  . 

Since  the  side  which  suffers  the  first  hypothetical  casualty  also 
suffers  the  first  actual  casuality  (at,  as  noted  above,  the  time  of 
the  first  hypothetical  casualty)  and  since  by  independence  and  non¬ 
atomicity  of  the  exponential  distribution 

p(ijj)jthe  T*(h,£)  and  T*(£,k)  are  all  distinct  from  one 
<27>  another}  =  1, 

it  follows  that 

(28)  P(i,j)|(Bv+,  Ry+)  e  j(i,j  -  1),  (i  -  l,j)f}  =  1 

and,  moreover,  that 

p(i,;i)l<Bv+>  Rv+)  =  (i’j  ■ 

_  p(i> j )|x*(k,X)  =  V*  for  some  k  and  £} 


i  j 

=  E  E 
k=l  1=1 


pC1? j )|T*(k, i)  =  V*} 


(we  have  used  (27)  here) 
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=  ij  P(l,j){T*(l,l)  =  V*f 

by  identical  distributions  of  the  T*(k,£). 

B 

Finally,  we  compute  that 
){T*C 1, 1)  =  V*} 

=  j  ^|t*(1,  1)  <  T*U,k)  for  k  <  i,  i  <  j  ; 

TJ(1,1)  <  T*(k,je)  for  2  <  k  <  i,  2  <  l  <  j; 

T*(l,l)  <  T*(l,i)  for  2  <  l  <  j; 

T*(l,l)  <  T*(k,  1)  for  .2  <  k  <  i} 

kl 

rijk2  +  (i  -  l)(j  -  l)kx  +  (j  -  l)kx  +  (i  -  l)kxJ  +  kL 

_ _ *1 _ 

ijk2  +  k^l  (i  -  l)(j  -  1)  +  (j  -  1)  +  (i  -  l)  +  1J 

=  kl 
ij(kx  +  k2) 

Thus 

p(1,j)«Bv+>  ^ =  13  k2y 


and  by  (28)  the  complementary  event  j(Bv+,  Ry+)  =  (i  -  J has 
5 ^-probability  given  as  follows: 
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L2  Homogeneous  Linear  Law  Process  With  Engagements 


To  study  this  process  we  must  introduce  further  notation, 
which  will  be  used  only  here.  Let 

E  =  N  x  N  x  N  , 

the  state  space  of  this  process,  be  endowed  with  the  discrete 
CT-algebra  £  and  the  discrete  topology.  We  take  as  the  sample 
space  n  the  set  of  mappings  m  =  (u^,  uj2  uj3)  of  [0,  ®)  into  the 
state  space  E  which  are  right  continuous  and  possess  left-hand 
limits  everywhere  (with  respect  to  the  discrete  topology  on  E  ). 

For  t  >  0  define 


B*(uu)  =  UU^t) 

R*(uO  =  iu2(t) 

and 

Dt(uu)  =  uu3(t)  . 

We  begin  discussion  of  the  process  with  binary  engagements  with  a 
model  of  a  three-component  process  in  which  B*  denotes  the  number 
of  unengaged  Blue  weapons  surviving  at  time  t,  R*  the  number  of  Red 
weapons  surviving  and  unengaged  at  time  t  and  the  number  of 
binary  engagements  (duels)  going  on  at  time  t.  Since  only  binary 
engagements  are  permitted,  Dt  is  then  also  the  number  of  surviving 
and  engaged  Blue  weapons  at  time  t,  which  necessarily  is  also  the 
number  of  Red  weapons  surviving  at  time  t  but  engaged  at  that 
instant. 

The  history  of  the  process  up  until  time  t  is  again  given  by 

Et  =  a((Bj,  RJ,  Du):  0  <  u  <  t) 

and  the  entire  history  by 

E  =  a(( BJ,  RJ,  Du):  u  >  0)  . 


107 


Let  3enote  tjie  probability  law  on  (fi,  F)  of  the  three 

component  process  ((B*,  R*,  D  ))  Q  governed  by  the  assumptions  of 
Process  L2,  subject  to  the  initiaT  conditions 

=  i  initially  unengaged  Blue  survivors 

Rg  =  j  initially  unengaged  Red  survivors 

Dq  =  k  engagements  initially  in  progress  . 

For  most  applications  one  would  have  k  =  0. 

(29)  THEOREM.  Under  the  family  of  assumptions  for  Process  L2, 
the  collection 

(0,  E,  Et,  (B*,  R*,  Dt),  P(l,;i,k)) 


is  a  regular  step  process  with 

a)  state  space  E  =  ; 

b)  jump  function  \  given  by 

Ui.j.k)  =  ku  +  ij (s.  +  s„)  ; 


c)  transition  kernel 

P((i,j,k),  (i  +  1,  j, 
P((i,j,k),  (i,  9+1, 
P((i, j ,k),  (i  +  1,  j  +  lr 
P((i, j ,k),  (i  -  1,  j  -  1, 


P  given  by 

Rp2u 

k  ~  ku  +  ij(s1  +  s2) 

kplU 

k  "  ku  +  ij(s1  +  s 2) 

kp3u 

k  ~  “  ku  +  ij(s1  +  s2) 

ij(sx  +  s2) 

k  +  ”  ku  +  ij(s1  +  s2)  ’ 
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d)  infinitesimal  generator  Q  given  by 


Q((i, j,k), 

Q((i,j,k), 


(i  +  l,  3,  k  -  1))  =  kp2u 
(  j  +  1,  k  -  1))  =  kpxu 
(i  +  1,  j  +  1,  k  -  1))-=  kp3u 


(i  -  1,  j  -  1,  k  +  1))  =  i j ( s x  +  s2) 
(i,j,h))  =  -  [ku  +  ijCs^  s2)]  . 


PROOF.  Denote  by  the  time  at  which  a  duel  is  first  initiated 
and  let  T2  be  the  time  at  which  the  first  initially  ongoing 
duel  terminates  (  =  +  ®  if  there  is  no  such  duel).  Then  the 
first  time  at  which  the  process  changes  state  is  given  by 


V  =  min{T15  T2  f  . 


From  the  proof  of  Theorem  (26)  and  assumption  7  we  see  that  with 
respect  to  the  probability  and  T2  are  independent 

and  exponentially  distributed  with  means  l/ij(s^  +  s2)  and  1/ku, 
respectively.  Therefore 

p(i,j,k)|v  >  >  t,  T2  >  tj 

=  P(i’j,k)jT1  >  t}P(i’j’k)|T2  >  t} 

=  exp  [  -  ij (s1  +  s2)t]  exp  (  -(ku)t] 

=  exp  [  -  \(i,  j ,k)t]  . 

At  time  V  exactly  one  of  four  events  occurs: 
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1)  One  engagement  terminates  with  the  .destruction  of  a  Red  weapon; 

2)  One  engagement  terminates  with  the  destruction  of  a  Blue 
weapon; 

3)  One  engagement  terminates  with  the  destruction  of  neither 
combatant; 

4)  One  new  engagement  is  initiated. 

If  the  current  state  is  (i,j,k)  then  the  new  states  reached  in 
these  four  cases  are 

1)  (i  +  1,  j,  k  -  1) 

2)  (i,  j  +  1,  k  -  1) 

3)  (i  +  1,  j  +  1,  k  -  1) 

4)  (i  -  1,  j  -  1,  k  +  1)  , 

respectively. 

By  now  familiar  computations, 


p(i,j>k)|(B*+)  r*+>  D*+)  =  (i  +  l,  j,  k  -  1)| 

=  <  ^  outcome  of  terminated  duel  is 

destroyed  Red  weapon} 

=  P(l»3»k){T2  <  TjJ  x 

p(i> j >k) joutcome  of  terminated  duel  is  destroyed 
Red  weapon} 


_ ku _ 

ku  +  IJTsT  +  s2) 


P 


2  * 


Similarly, 

RJ+>  Dv+)  =  (i,  j  +  1,  k  -  1)| 
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_  p(i>J>k)j^  <  T^,  outcome  of  terminated  duel  is  destroyed 
Blue  weapon } 

=  P^1,:^k^{T2  <  Tx}  x 

p(i> j {outCOme  is  destroyed  Blue  weapon} 

_ ku _ 

~  ku  +  ij(s1  +  s  2 )  P1  » 


and 


pCi,  j,K)  {(B*+>  R*+,  Dy+)  =  (i  +  1,  j  +  1,  k  -  1)) 

=  <  T^,  outcome  of  terminated  duel  is  destruction 

of  neither  weapon  } 

=  P(l,j,k)  fr2  <  Tx}  x 

p( j »k)  {outcome  is  mutual  survival] 


ku 


ku  +  ij(s1  +  s2)  p3 


P7  • 


Finally, 

pCi,j,k){(B*+>  r*+>  Dy+)  =  (i  -  1,  j  -  1,  k  +  1)} 
_  p(i,  j  (T^  <  T^} 

N 

i-Ks-t  +  s2) 

-  ku  +  IJTsT  +  s2) 

Standard  arguments  now  complete  the  proof. 
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Consider  now  the  two-dimensional  process  ((Bt,  Rt))t>0  of  total 
numbers  of  survivors  defined  by 


Bt  =  B*  + 


and 


Rt  =  R*  + 


Unfortunately,  ((Bt,  Rfc))  q  is  not  a  Markov  process  with  respect  to 
the  u-algebras  (Et)t>Q  because  for  each  t  the  history  F  of  the 
process  up  until  time’  t  contains  information  about  the  process, 
namely  the  components  B*,  R*  an<^  which  is  not  contained  in  the 
pair  (Bt,  Rt). 

Functions  of  Markov  process  have  been  studied  but  mostly  in  terms 
of  limiting  behavior  (e.g.,  Laws  of  Large  Numbers  and  Central  Limit 
Theorems,  cf.  Chung  (1967))  and  identifiability  problems  (cf. 
Rosenblatt  (1971))and  hardly  at  all  in  terms  of  transient  behavior 
and  computation  of  distributions.  Further  research  is  desirable. 

We  can  make  some  simple  remarks  about  computations.  If  for  each 
t  we  know 

Pt((i,j,R),  U,m,n))  =  P(lj:)’k)|(B*,  R*,  Dt)  =  (£,m,n)[ 
for  all  (i,  j,  k)  and  (£,  m,  n),  then  we  can  compute 

=  E  Z  Pt((iJ>R)5  (i,m,n))  . 

m  n+£  =  |u 


This  computation  is  not  so  hard  as  appears  at  first  glance,  since  we 
know  from  the  theory  of  regular  step  processes  that 


Pt((i,j,k),  (£,m,n))  =[exp  (tQ)]  ((i,j,k),  (£,m,n)), 
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and  hence  can  be  approximated  as 


M  ^p 

P^-CCij j jk)5  (A,m,n))  ~  E  — r  Qp((i,j,k)>  (A,m,n)) 

r  p=0  p* 

Since  by  linearity  of  expectation  we  have  for  each  (i,j,k)  e  E 
(-3^]  =  E(i>j,k)  [B*]  +  [Dt] 

any  scheme  (in  particular,  that  discussed  in  Section  IV  of  the  paper) 
for  approximating  expectations  of  a  regular  step  process  also  approxi¬ 
mates  expected  numbers  of  survivors  in  this  process. 

To  produce  a  good  probabilistic  characterization  of  the  process 
((Bt,  R  ))  will  take  further  research  into  the  behavior  and  charac¬ 
terization  of  functions  of  Markov  processes. 


L3  Heterogeneous  Linear  Law  Process 


Assume  that  there  are  M  Blue  weapon  types  and  N  Red  weapon 
types.  The  state  space  E  for  the  heterogeneous  linear  law  process 
is  then  given  by 

E  =  N  x  ...  XN(M+N  times)  , 

along  with  the  discrete  a-algebra  |  .  As  sample  space  Cl  we  have 
the  family  of  functions  from  [0,  ®)  into  E  which  are  right- 
continuous  and  have  left-hand  limits  everywhere.  For  w  =  (to-^, 

^M’  ^1’  ’  *  *  ’  w^ere  Wp  ...j  uu^  are  mappings  of 

[0,  oo )  into  which  are  right-continuous  and  have  left-hand  limits, 
with  respect  to  the  discrete  topology  on  N  ,  and  t  >  0,  we  define 

Bt(u))  =  (u>j(t),  ...,  t:) ) 

Rt(aj)  =  (ou£(t),  ...,  oujJ(t)) 

B^(uu)  =  to^t)  ,  i  =  1,  ...,  M  , 

R^(w)  =  uuj(t)  ,  j  =  1,  ...,  N  . 

is  the  number  of  type  i  Blue  weapons  surviving  at  time  t,  Bt  = 

M  i 

(B',  ...,  B  )  the  Blue  force  surviving  at  time  t,  R"  the  number  of 

L  T- 

type  j  Red  weapons  surviving  at  time  t  and  R  the  surviving  Red  force 
at  time  t. 

The  history  of  the  attrition  process  up  until  time  t  is 
It  =  a((§s’  ~s):  0  <  s  < 

and  the  entire  history  is 

E  =  crC(Bs,  Rs):  s  >  0)  . 
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Points  in  E  are  written  in  the  form  (x,y)  with  x  e  NM  and  x. 

N  ~  i 

denoting  the  Blue  type  i  weapons  and  y  e  N  with  y .  the  number  of 
type  j  Red  weapons . 

For  each  (x,y)  e  E  let  P^X,y^  be  the  probability  law  on  (fi,F) 
of  the  attrition  process  governed  by  the  family  of  hypotheses  of 
Process  L3,  subject  to  the  initial  conditions 


and 


x 


(30)  THEOREM.  Under  the  assumptions  of  Process  L2,  the  process 

(Cl,  E,  Et,  (Bt,  Rt)  ,  P(x’y)) 


is  a  regular  step  process  with 

a)  state  space  (E,  E)  as  defined  above; 

b)  jump  function  \  given  by 

M  N 

(31)  \(x,y)  =  E  £  x±y  [k1(j,i)  +  k2(i,j)] 

i=l  j=l  3 

where  k  (q,v)  =  s  (u,v)p  (n,v)  for  q  =  1,  2  and  appropriate  values 
of  |j  and  v  ; 

c)  transition  kernel  P  given  by 


P((x,y);  (x;  yv 


yj-1’ 


yN^ 


M 

y  E  k  (j,i)x. 
3  i=l  1  1 

l(x,y) 


(32) 

P((x,y);  ( x^j 


x^  I,  •  •  •  j  x^;  y)) 


N 

x.  E  k?(i,j)y. 
3=1  3 

X(x,y) 
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d)  infinitesimal  generator  Q  given  by 


M 


Q.C  ( x,  y ) ,  (xj  Yjj  •••5  Yj  ~  ))  —  Yj  ^  k^(j,i)x^ 


M  N 


Q((x,y),  (x,  y) )  =  -  E  E  x.y.(k  Q,i)  +  k„(i,j 

i=l  j=l  1  0  x  * 


)) 


N 

Q((x,y),  (Xp  ...,  x±  -1,  ...,  xM;  y))  =  x±  k2(i,j)y..  . 

PROOF.  We  follow  the  pattern  established  in  the  proof  of  Theorem  (26), 
only  now  the  notation  is  more  complicated. 

+“A_ 

Denote  by  T*(i,j;  k,£)  the  time  at  which  the  k^  Blue  weapon  of 
•D  til 

type  i '.first  hypothetically  detects  and  kills  the  Red  weapon  of 

type  j  and  by  T*(j,i;  £,k)  the  time  at  which  the  first  hypothetical 

kill  of  that  specific  Blue  weapon  by  the  same  particular  Red  weapon 

occurs.  With  respect  to  the  probability  law  P^X’^,  for  each  i 

•k 

and  j  the  random  variables  T|(i,j;  1,1),  Tg(i,j;  1,2),  ..., 

T*(i,jj  l,yj),  T*(i,jj  2,1),  ...,  T*(i, j ;  2,^),  ...,  T*(i,j;  xi,l), 
...,  Tg(i, j ;  x^,y^)  are  independent  and  identically  exponentially 
distributed  with  mean  l/k1(j,i).  Similarly,  T*(j,i;  1,1),  ...,  T* 

(j,i;  l,xi),  ...,  T*(j,i;  yj}l),  ...,  T*(j,i;  yj?  x±)  are  independent 

and  identically  exponentially  distributed  with  mean  l/k9(i,j).  These 

( x  vO 

two  families  of  random  variables  are  independent  with  respect  to  Pv  ’ J ' 
and  as  the  weapon  types  i  and  j  vary  all  families  of  random  variables 
so  obtained  are  independent. 

With  respect  to  p(x’^  the  time 


V*  =  min  min  min{T*(i,  j  ;  k,£);  T*(j,i;  £,k)| 

i=l, . . .  ,M  k<x . 

j=l, . ..,N  "  1 

£<yj 
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of  the  first  hypothetical  fatality  and  the  time 


V  =  inf {t :  (Bt,  Rt)  t  (BQ,  RQ)} 


at  which  the  first  actual  fatality  occurs  are  equal  almost  surely. 
By  the  independence  hypothesis  5,  we  have 

p(x,y)j-v*  >  t} 

_  p(x> y'{T*(i, j ;  k,£)  >  t,  T*(j,i;  £,k)  >  t  for  all 

appropriate  i,  j,  k, 


=  n  n  p^x’y){T*(i,  j ;  k,£)  >  t  for  all  k  <  x. ,  JL  <  y -3 
i=l  j=l  J 


x  n  n  p(x>y){T*(j,i;  £,k)  >  t  for  k  <  x. ,  l  <  y .} 
j=l  i=l  K  J 

M  N  Xi  yj  fxy. 

=  n  n  n  n  ,y;[T*(i,j;  k,*)  >  t} 

i=l  j=l  k=l  SL=1 


N  M  1  1  (y 

x  n  n  n  n  pc  ,y'  {T*(j,i;  z,k)  >  t} 

j=l  i=l  1=1  k=l 


=  n  n  n  n  exp  [  -  k1(j,i)t] 
i=l  j=l  k=l  1=1 


N  M  Jj  l 

x  n  n  n  n  exp  [  -  k2(i,j)t] 

j=l  i=l  1=1  k=l 


M  N 

=  n  n  exp  [  -  x.y.k  (j,i)t] 
i=l  j=l  J 


N  M 

x  n  n  exp  [  -  y.XjkgCijjJt] 
j=l  i=l  J 
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M  N 

=  exp  [  -  t  z  Z  x.y.(k  (j,i)  +  k«(i,j))] 
i=l  j=l  1  3  1  * 

=  exp  [  -  X(x,y)t  ]  , 
where  \  is  defined  by  (31),  so  that 

P(X’y){V  >  t}  =  exp  [  -  X(x,y)t] 

as  required. 

Since  almost  surely  with  respect  to  P^x,y^  the  T*(i,j;  k,£)  and 

B 

T^(j,i;  £,k)  are  all  distinct  from  one  another,  it  follows  first  of 
all  that 

P  5  {(By_j_,  ^  yi-"^’  y2*  ’  ’  ’  ’  )  5  •  ••> 

(xj  yv  •  ••>  yjj-i)j  ^^i  —  •  •  •>  j  y), 

. . . ,  ( x^ ,  . . . ,  x^~ ^"5  y ) } }  —  i  • 

For  a  fixed  j  we  have 


p(  —  (x>  yv  •  yj~^->  ***>  yN)i 


=  p^x>^){v*  =  T*(i,j;  k,£)  for  some  i,  k  <  x.;  and  &  <  y.} 
B  —  l  J 


x/  x .  y . 
M  2-  3 


=  Z  2  2~  P(X,y){V*  =  T*(i,j;  k,£)} 

i=l  k=l  1=1  1  -  B 

=  y.  2  x.  p(x>y){v*  =  T*(i,j;  1,1)}  • 


3  i=l  1 


Bv 


But 


p(x,y)fv*  =  1,1)3 
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=  P(X,y){T*(i,j;  1,1)  <  T*(j,i;  £,k)  for  all  k,  ly 

T*(i,j;  1,1)  <  T*(n,m;  £,k)  for  m  /  i,  n  ? j ,  all  k,  £ 

T*(i,j;  1,1)  <  T*(m,n;  k,£)  for  m  /  i,  n  1  j,  all  k,£ 

T*(i,j;  1,1)  <  T*j(i,j;  p,q)  for  p  >  1  or  q  >  l) 

^y^kgCijj)  +  E  E  x  y  [k2(m,n)  +  k  (n,m)]  +  x.y.k  (j,i) 

J  m^i  n^j  J 

=  P((x,y);  (x;  y±,  ...,  y^-1,  ...,  yN)) 

where  P( • ,  • )  is  defined  by  (32).  The  proof  that 

p(  Ry_^_)  —  (X^,  •••,  X^~l,  •••>  y)} 

=  P((x,y),  (xx,  xi-l,  xM;  y) 

for  each  i  is  entirely  analogous  and  is  hence  omitted. 

The  proof  is  complete  upon  application  of  standard  arguments 
based  on  the  ’’memoryless”  property  of  the  exponential  distribution 
which  extend  the  preceding  computations  to  show  that  if  one  defines 
for  each  t  the  time 

Vt  =  inf {s  >  t:  (Bs,  Rg)  /  (Bt,  Bt)} 

of  the  first  transition  of  the  process  after  time  t  then 
P(x’y){Vt  -  t  >  u  |  Et]  =  exp  [  -  X(Bt,  Rt)u] 

and 

P(X’y)t(BVt+,  R v  )  =  a)  =  P((Bt,  Rt);  a) 
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for  all  (x,y)  and  a  e  E.  From  these  two  statements  one  concludes 
that  (n,  E,  Et>  (Bt,  gt),  is  indeed  a  regular  step  process 

with  jump  function  X  and  transition  kernel  P. 


Ml  Hetergeneous  Mixed  Law  Process 


For  simplicity  and  clarity  we  assume  that  each  side  possesses 
exactly  two  types  of  weapons,  one  of  which  obeys  assumptions  of 
weapons  in  square  law  models  and  the  other  of  which  behaves  accord¬ 
ing  to  the  assumptions  of  the  linear  law  models.  Our  convention  is 
that  Type  1  weapons  satisfy  square  law  assumptions  and  Type  2  weapons 
linear  law  assumptions. 

4 

The  state  space  becomes  E  =  N  with,  as  usual,  the  discrete 
a-algebra  E  .  As  sample  space  Q  we  have  the  family  of 
right-continuous  functions  from  [0,  ®)  into  E  which  everywhere 
have  left-hand  limits.  For  w  =  (u)p  <*>4)  e  0  we  define 

B^(uj)  =  u^Ct) 

B^(w)  =  «u2(t) 

rJXw)  =  u)3(t) 

R^(u0  =  uu4(t)  . 

For  q  =  1,  2,  (resp.,  R^)  is  the  number  of  type  q  Blue  (resp., 

Red)  weapons  surviving  at  time  t.  We  further  define 

St  =  <®t>  Bt> 

and 

St  =  (Rt> 

which  are  the  surviving  Blue  and  Red  forces,  respectively,  at 
time  t. 

As  usual  the  history  of  the  process  until  time  t  is 
Et  =  c((Bg,  Rs):  0  <  s  <  t) 
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and  the  entire  history  is 


E  =  a((Bs,  Rg):  s  >  0)  . 

( i  i  k  £) 

Denote  by  Pv  ’  ’  ’  '  the  probability  law  on  the  measurable 
space  (Q,E)  of  the  attrition  process  governed  by  the  family  of 
eight  assumptions  of  the  Process  Ml,  conditioned  on  the  event 

{B0  =  i’  B0  =  j’  R0  =  k>  Ro  =  * 

(33)  THEOREM.  Under  the  assumptions  of  Process  Ml, 

(0,  E,  Et,  (Bt,  Bt),  P(i,j,k,'e)) 


is  a  regular  step  process  with 

a)  state  space  (E,  E),  as  defined  above; 

b)  jump  function  X  given  for  nonabsorbing  states  [that  is, 
states  (m,n,p,q)  such  that  m  +  n  >  0  and  p  +  q  >  0]  by 

^'kp1(  1)  +  AP1(2)^ 


X(i,j,k,je)=  ir^ - ^ 


+  j(ks1(l)ql(l)  +  £s1(2)q1(2)) 
^iP2d)  +  JP2C2)  ^ 

+  Rr2  \ — m - j 


+  £(is2(l)q2(l)  +  js2(2)q2(2))  ; 

c)  transition  kernel  P  given  by 


P((i» j »k,£);  (i, j ,  k-1, £))  = 


P((i,j,k,A);  (i, j ;  k,4-l)) 


iripi(1)  wt +  jk  *  s1(1)q1(1) 

X(i,j,k,je) 

irlPl(2)  I+lc  + 

X(i,j,k,X) 
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? 


kr9p2(l)  —  +  li  •  s2(l)q2(l) 

p((i,j,k,A);  (i-i,j,k,D)  = — — — itt - £ — 


P((i,j,k,£);  (i,j-l,k,£)) 


4r2P2(2)  -jJj  +  £j  •  s2(2)q2 

\'(i,  j,k,X  ; 


(2) 


d)  infinitesimal  generator  Q  given  for  nonabsorbing  states 
(i,j,k,j&)  by 


•Q(  (i>  j  X)>  (ij 3  ,k-l,A))  -  P-^Cl)  +  jk  *  s^(l)q^(l) 

.Q((ij j jk,£);  (i,j,k,X-l))  =  irr  P1(2)  ^  +  U  *  s1(2)q1(2) 

(  /kp  (1)  +  £p,(2)\ 

Q((i,j,k,x);  (i,j,k,i))  =  -|ir1^  ■■  -y-— - j 

+  j(ks1(l)q1(l)  +  £s1(2)q1(2)) 
/ip2(l)  +  jp2(2)\ 

+  kr2  { — m — ) 

+  £(is2(l)q2(l)  +  js2(2)q2(2)) 
Q((i, j,k, X);  (i-l,j,k,X))  =  kr2p2(l)  ^  +  Xi  *  s2(l)q2(l) 

Q((i>j>k,X);  (i, j-l,k,jO)  =  kr2p2(2)  -j^*-  +  Xj  •  s2(2)q2(2)  . 

If  state  (i,j,k,x)  is  absorbing  then Q((i,j,k,X),  a)  =  0  for  all 
a  e  E. 

PROOF.  For  one  last  time  we  resort  to  the  fiction  of  hypothetical 
kills. 

Denote  by  T*(m,n)  the  time  at  which  ■some  Red  type  n  weapon  is 
B 

first  hypothetically  killed  by  some  Blue  type  m  weapon  and  by 
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T*(n,m)  the  time  at  which  some  Blue  type  m  weapon  is  first 
hypothetically  killed  by  some  type  n  Red  weapon,  for  m,  n  =  1,  2. 


Hypothetical  shots  fired  by  a  type  1  Blue  weapon  form  a  Poisson 
process  with  rate  rp,  so  such  shots  directed  at  type  1  and  type  2 
Red  weapons  form  independent  Poisson  processes  with  rates 


k  <L 

and  r^  ,  respectively,  relative,  of  course,  to  the 


f  i  i  k  £ ) 

probability  measure  Pv  ,J>  ’  .  Hypothetical  fatal  shots  fired  by 

this  particular  Blue  type  1  weapon  against  type  1  and  type  2  Red 
weapons  thus  form  independent  Poisson  processes  with  respective 


k  SL 

rates  r^p^(l)  and  r^p^(2)  and  hence  by  assumption  8  and 
the  Superposition  Theorem  for  Poisson  processes  it  follows  that 


p(i,J,k,£)^T*(ljl)  >  =  exp[  _  t(irlPl(l)  ^)] 

(34) 

P(l,:’k,£){T|(1,2)  >  t}  =  exp  [  -  t(irlPl(2)  )]. 

Entirely  analogously  we  have 

p(i,j,k,£)£T*(ljl)  >  =  exp  [  •_  t(kr2p2(l)  jij)] 

(35) 

p(1, 3 > £ )[T*(1} 2 )  >  t]  =  exp  [  -  t(kr2p2(2)  Uj)]  . 

Hypothetical  detections  of  a  particular  type  n  Red  weapon 
by  a  particular  type  2  Blue  weapon  form  a  Poisson  process  with 
rate  s^n),  and  hypothetical  fatal  shots  by  the  same  Blue  weapon 
at  the  same  Red  weapon  a  Poisson  process  with  rate  s1(n)qp(n)  from 
which  one  obtains  by  the  independence  assumption  8  that 

P(i,j,k,£){T*(2,1)  >  t}  =  exp [  -  t(jks1(l)q1(l)] 

(36) 

P(i,j,k,£){T*(2,2)  >  t}  =  exp [ . -  t(j£s1(2)ql(2)]  . 
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Similarly 


l  *? 


p(i,j,k,£)|T*(2jl)  >  tJ  =  exp[  _  t(£i  s2(l)q2(l))] 

(37) 

p(i,3,k,i)jT*(2j2)  >  t|  =  exp  [  -  t(£j  s2(2)q2(2))]. 

But  the  independence  assumption  8  implies  that  with  respect  to 

p(i>j>k,£)  random  variables  T*(m,n),  m,  n  =  1,  2  and  T*(n,m), 

•D  R 

n,  m  =  1,  2  are  mutually  independent.  The  minimum  V*  of  these  eight 
random  variables  is,  moreover,  equal  with  probability  one  to  the 
time 


V  =  inf j t:  (Bt,  Rt)  +  (BQ,  RQ) }  . 


Therefore,  by  (34),  (35),  (36),  and  (37) 

2  2.. 

p(i,j,k,£)|v  >  t}  =  n  n  p(157}k,£r)|T*(m,n)  >  t}p(  >k>£^T*(n,m)  >  t 

m=l  n=l 

=  exp  t  -  t(irlPl(l)  ^  ♦  ir2pl(2) 

*  kr2P2(1)I?J+  kr2p2(2)iij 

+  jk  •  s1(l)q1(l)  +  ji  •  s1(2)ql(2) 

+  ii  •  s2(l)q2(l)  +  £j  •  S2(2)q2(2))] 

=  exp  [  -  tX(i,j,k,£)] 

by  a  simple  rearrangement  of  terms. 

With  probability  one  the  eight  aforementioned  random  variables  are 
all  distinct  from  one  another.  We  then  have,  for  example, 
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p(i5j,k5£)l(BV+J  R^)  = 

=  P(l,:i,k’£)|v*  =  T*(  1, 1)  or  V*  =  T*( 2, 1)| 

=  P(l,:i,k’A)iV*  =  T*(l,l)}  +  jv*  =  T*(2,i)| 

=  P(l’j>k>X){T*(l}l)  <  min{ T*(  1, 2),  T*(2,l)  , 

T*(2,2),  T*(l,l)  , 

Tr(  1’  2)>  T*(2,l)  , 

Tr( 2  j  2 ) }  f 


+  P 


|T*C 2 , 1)  <  minjTj(l,l),  T*(l,2)  , 


Tg( 2 j  2  ) ,  T*( 1,1)  , 
Tr(  1}  2)j  T*(2,l)  , 


Tg(2,2)| 


irlpl(1)  k+£  jk  •  s.CDq.Cl) 

-  +  - - - 


x(i,  j,k3je) 


\(i, j,k,£) 


here  we  have  omitted  some  long  and  unenlightening  expressions  occurring 
in  the  course  of  the  computation. 


One  similarly  checks  that 

P(i’j,k’‘)|(Sv+,  ?v+)  =  (i,j,k,£-l)! 

lrlpl(2)  +  '  S1<2)<:|1<2> 


l(ijj  jkj  jj ) 
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that 


P(1’j’M>|CB,v+,  Rv+)  =  (i-1,  j,  k,  i)| 


kr2p2(l)  +  iJL  .  s2(l)q2(l) 
X(i, j,k,£) 


and,  finally,  that 

P(i,j,k,£)|(Bv+,  Rv+)  =  (i,j-l,k,A)f 

kr2p2C2)  jjj  +  j4  *  s2(2)q2(2) 

Ui,  j,k,£) 

Hence  in  all  cases  we  have 

P(l,3,k,£){(Bv+>  Rv+)  =  a\=  P((i,j,k,jO;  a) 


and  the  proof  of  the  Theorem  is  complete  following  standard  exten¬ 
sions  based  on  the  "memory less"  property  of  the  exponential 
distribution. 


Mia  Homogeneous  Mixed  Law  Process 


We  assume  that  the  Blue  side  has  only  type  1  weapons  and  the 
Red  side  only  type  2  weapons. 

Consider  the  4-dimensional  Heterogeneous  Mixed  Law  Process  Ml 

derived  above.  With  respect  to  a  probability  of  the  form 
2  1 

and  are  identically  zero  almost  surely  and  can  be  disregarded. 

Now 


X(i,0,0,£)  =  ir1P1(2)  +  i£  .  s2(l)q2(l) 


while 

ir  p  (2) 

P((i,0,0,i),  (i,0,0,i-l»  x(i}o,0,j£) 


and 

i£  •  s2(l)q2(l) 

P((i,0,0,£),  (i-l,0,0,£))  =  x(i,0,0,£) 


Hence  we  have  the  following  result. 


(38)  THEOREM.'  Let  the  assumptions  of  Process  Ml  be  satisfied. 
Then  the  stochastic  process 


(0,  E,  Et,  (Bj,  R2t), 


is  equivalent--in  the  sense  of  being  a  Markov  process  with  the  same 
infinitesimal  generator  (and  thus  a  regular  step  process  with  the 
same  jump  function  and  transition  kernel) — to  any  regular  step 

A 

process  with  state  space  N  x  N,  jump  function  X  given  by 
X(i,£)  =  ir1p1(2)  +  i£  •  s2(l)q2(l) 
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and  transition  kernel  P  given  by 


ir  p  (2) 

±i  .  s9(l)q9(l) 

P((i,£);  (i-l,£))  =  x(i,£)  ' 

From  Theorem  (38)  one  trivially  obtains  our  final  Theorem. 


(39)  THEOREM.  The  attrition  process  governed  by  the  family  of 
assumptions  of  Process  Mia  is  a  regular  step  process  with 

a)  state  space  N  x  N; 

b)  jump  function  \  given  by 


X(i,j  ) 


irp  +  ijsq 

0 


if  i  >  0,  -j  >  0 
if  i  =  0  or  j  =  0  ; 


c)  transition  kernel  P  given  for  states  (i,j)  with  i  >  0  and 
j  >  0  by 


P((i,j),  (i,j-l)) 


rp 

rp  +  jsq 


P((i,j),  (i,j-l)) 


jsq  . 
rp  +  jsq  ’ 


d)  infinitesimal  generator  Q  given  for  states  (i,j)  with  i  >  0 
and  j  >  0  by 


Q((i,j),  (i,j-l))  =  irp 
Q((i,j),  (i> j  ))  =  -  (irp  +  ijsq) 

Q((i,j),  (i-1,  j  )  )  =  ijsq  . 

If  i  =  0  or  j  =  0,  Q((i,j);  (k,£))  =  0  for  all  (k,£)  e  E.  [j 


The  particularly  simple  form  of  the  transition  kernel  P  is 
noteworthy.  Each  column  is  constant;  perhaps  this  structure  can  be 
exploited  in  computational  applications . 
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